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developed  to  investigate  the  steady  state  performance  characteristics  of  the  heat 
exchanger  over  a  wide  range  of  air  inlet  temperatures,  mass  flow  rates,  and  reser¬ 
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cycle.  Computational  results  are  shown,  and  a  brief  comparison  with  the  experiments 
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ABSTRACT 


This  report  documents  a  study  of  a  prototype  air- 
to-air  heat  exchanger  of  which  the  thermal  conductivity  is 
governed  by  the  reservoir  temperature  of  the  heat  pipes.  The 
work  is  introduced  with  a  brief  review  of  the  basic  operating 
principles  of  heat  pipes  and  of  the  various  existing  schemes 
for  variable  conductance. 

An  experimental  facility,  incorporating  a  compu¬ 
terized  data  acquisition  system,  was  developed  to  investigate 
the  steady  state  performance  characteristics  of  the  heat 
exchanger  over  a  wide  range  of  air  inlet  temperatures,  mass 
flow  rates,  and  reservoir  temperatures.  Some  tests  were  con¬ 
ducted  imposing  a  single  operating  temperature  on  all  of  the 
heat  pipes,  while  in  others  the  heat  pipe  rows  were  allowed  to 
operate  independently.  The  results  are  expressed  in  terms  of 
heat  exchanger  effectiveness  versus  normalized  reservoir  tem¬ 
perature  for  various  values  of  the  ratio  of  the  mass  flow 
rates  of  the  air  streams.  Conclusions  are  drawn  as  to  the 
internal  behaviour  of  this  type  of  heat  pipes  under  various 
operating  regimes. 

An  original  computer  model  is  presented  which  com¬ 
bines  a  quasi-analytical  solution  for  the  air  and  fin  tempera¬ 
ture  profiles  in  the  heat  exchanger  with  a  heat  pipe  simulation 
routine  based  on  the  hydrodynamic  equations  involved  in  the 
working  fluid  cycle.  Computational  results  are  shown,  and  a 
brief  comparison  with  the  experiments  is  made.  Possible 
modifications  that  would  imporve  the  accuracy  and  the  versi- 
tility  of  the  program  are  discussed. 
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1.  INTRODUCTION 


A  problem  that  has  always  been  associated  with  the 
design  of  effective  air-to-air  heat  exchangers  is  their  high 
bulk  and  weight,  due  to  the  very  large  surfaces  required  to 
transfer  significant  amounts  of  heat  from  one  airstream  to 
another.  The  birth  and  rapid  growth  of  heat  pipe  technology 
in  the  1960s  led  to  the  eventual  development  of  more  compact 
heat  exchangers  taking  advantage  of  the  very  high  thermal  con¬ 
ductivity  of  these  devices.  The  continued  development  of  such 
heat  exchangers  brought  the  promise  of  many  new  applications 
in  the  fields  of  energy  conservation  and  air-conditioning.  At 
the  same  time,  however,  a  significant  limitation  of  existing 
heat  pipe  heat  exchangers  became  apparent,  namely,  a  lack  of 
controllability.  While  it  is  undisputably  a  high  performance 
heat  transfer  device,  the  standard  heat  pipe's  thermal  con¬ 
ductivity  is  highly  dependent  on  its  operating  temperature  and 
on  the  thermal  load  placed  across  it;  in  other  words,  the  over¬ 
all  heat  transfer  coefficient  Is  a  function  of  the  heat  sink 
and  source  conditions.  With  the  standard  heat  pipe,  it  is 
impossible  to  adjust  this  heat  transfer  coefficient,  and  there¬ 
fore  impossible  to  regulate  the  heat  flux  in  a  heat  exchanger 
using  them. 

This  limitation  poses  problems  for  heat  recovery  and 
temperature  control  applications.  In  heat  recovery  systems, 
in  which  waste  heat  is  recovered  from  exhausted  air  (or  process 
gases)  and  used  to  heat  cold  incoming  air,  it  is  desirable  to 
maximize  the  heat  transfer  rate,  i.e.,  the  heat  exchanger's 
effectiveness.  Uncontrollable  heat  pipe  heat  exchangers  work, 
however,  at  an  effectiveness  which  is  not  necessarily  their 
optimum.  The  working  point  is  passively  determined  by  the 
temperature  and  mass  flow  rate  of  the  two  airstreams.  Such  heat 
exchangers  are  also  inadequate  to  meet  the  requirements  of  air- 
conditioning  systems,  where  one  must  control  the  exit  tempera¬ 
ture  of  one  of  the  ai rstreams,  si  nee  the  air  exit  conditions  are 
uniquely  determined  by  the  inlet  conditions.  This  lack  of 
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controllability  is  the  motive  behind  the  research  on  variable 
conductance  heat  pipe  heat  exchangers.  Much  work  has  been 
done  to  develop  both  passively  and  actively  controlled  heat 
pipes.  Some  of  the  resulting  systems  will  be  mentioned  in 
the  following  sections. 

This  report  documents  both  an  experimental  and  nurneri 
cal  study  of  a  prototype  variable  conductance  heat  pipe  (VCHP) 
heat  exchanger  based  on  the  principle  of  operating  temperature 
modulation.  This  particular  design,  which  makes  use  of  liquid 
reservoirs  of  excess  working  fluid,  offers  the  possibility  of 
heat  flux  maximization  and  temperature  control  without  moving 
parts.  The  objectives  of  this  study  are  to  determine  the  steady 
state  performance  characteri sti cs  of  this  he_t  exchanger  over 
a  wide  range  of  operating  conditions  and  to  apply  heat  pipe 
theory  to  develop  a  numerical  model  of  the  heat  exchanger 
capable  of  predicting  temperature  profiles  and  heat  flux  for 
specified  air  inlet  conditions.  Such  a  model  would  be  useful 
as  a  tool  for  both  the  design  and  control  of  this  type  of  heat 
exchanger. 
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2.  HEAT  PIPE  PRINCIPLES 

2 . 1  The  standard  heat  pipe 

Before  discussing  the  details  of  the  liquid  reservoir 
VCHP,  it  is  useful  to  review  the  fundamental  principles  of  heat 
pipes.  The  heat  pipe  is  a  high  performance  heat  transfer  device, 
able  to  transfer  heat  at  high  rates  over  considerable  distances 
with  almost  no  temperature  drop.  In  its  basic  form  ti'.e  heat 
pipe  consists  of  a  sealed  tube  lined  with  a  porous  capillary 
wick  that  is  saturated  with  a  volatile  fluid  as  shown  in 
figure  1.  Heat  absorbed  by  the  evaporator  zone  vaporizes  the 
fluid;vapour  flows  to  the  other  end  where  it  condenses,  giving 
up  the  heat.  The  fluid  cycle  is  then  completed  as  the  condensed 
fluid  flows  in  the  wick  back  to  the  evaporator.  The  axial  heat 
transfer  is  thus  equal  to  the  mass  circulation  rate  times  the 
latent  heat  of  vaporization  of  the  working  fluid. 

Frequently  the  evaporator  and  condenser  sections  are 
separated  by  an  adiabatic  zone,  where  there  is  no  heat  flux 
into  or  out  of  the  pipe.  A  heat  pipe  can,  in  principle, 
operate  with  the  evaporator  above  the  condenser  as  shown  in 
figure  1,  in  which  case  the  capillary  head  generated  in  the 
wick  must  overcome  the  gravitational  head,  but  in  many  cases 
the  evaporator  is  intentionally  placed  below  the  condenser  so 
that  these  two  forces  work  together,  which  facilitates  the 
return  of  the  condensate  to  the  evaporator. 

What  appears  to  be  a  simple  device  actually  involves 
a  complex  set  of  thermal  and  hydrodynamic  processes,  which 
have  made  it  difficult  for  researchers  to  draw  firm  conclusions 
about  the  exact  nature  of  the  heat  transfer  processes  in  a 
particular  heat  pipe  for  a  given  set  of  heat  source  and  sink 
conditions.  What  has  been  determined  is  that  the  heat  pipe's 
performance  is  bounded  by  an  envelope  of  limiting  conditions 
as  shown  in  figure  2. 
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2 . 2  Limits  to  heat  transfer 

The  range  of  vapour  temperature  (also  called  opera¬ 
ting  temperature)  over  which  a  particular  heat  pipe  functions 
is  determined  principally  by  the  working  fluid  used.  At  the 
kv  and  of  this  range,  fluid  viscosity  imposes  a  limit  on  the 
heat  transfer  capacity  of  the  heat  pipes,  as  shown  in  figure  2. 

For  slightly  higher  operating  temperatures,  what  is 
called  the  sonic  limit  is  the  dominant  limiting  factor.  It 
can  be  shown  analytically  that  there  is  a  correspondence  be¬ 
tween  constant  area  flow  in  a  heat  pipe  with  mass  addition 
(evaporator)  and  removal  (condenser)  and  constant  mass  flow 
in  a  converging-diverging  nozzle.  Consequently,  just  as  there 
is  a  sonic  (Mach  =  1)  limit'tion  on  the  flow  velocity  through 
a  nozzle,  there  is  a  similar  limit  on  the  flow  velocity  through 
the  adiabatic  section  of  the  heat  pipe. 

Another  phenomenon  which  may  limit  heat  flux  over  a 
certain  range  of  operating  temperature  is  liquid  entrainment 
in  the  vapour.  In  a  heat  pipe,  the  vapour  and  the  liquid 
generally  flow  in  opposite  directions.  Since  they  are  in 
contact  at  the  wick  surface  (if  the  wick  is  saturated),  this 
sets  up  a  mutual  drag  at  the  vapour-liquid  interface.  If  the 
relative  velocity  between  the  liquid  and  vapour  becomes  too 
great,  the  interface  becomes  unstable  and  droplets  of  liquid 
will  be  torn  from  the  wick  and  entrained  in  the  vapour.  Since 
this  liquid  never  reaches  the  evaporator,  it  cannot  contribute 
to  the  heat  transferred  by  the  heat  pipe.  Thus,  the  maximum 
axial  heat  transfer  is  no  longer  equal  to  the  maximum  fluid 
circulation  rate  times  the  latent  heat  of  vaporization,  but 
is  some  lower  value. 

The  fact  that  there  exists  a  maximum  capillary  head 
for  any  wick-fluid  combination  results  in  a  hydrodynamic  limit 
on  heat  pipe  capacity.  In  order  for  the  heat  pipe  to  operate, 
the  following  pressure  balance  equation  must  be  satisfied 


between  any  two  points  along  the  heat  pipe  : 


‘Pc.raax  +  ‘Pg  >  4Pt  +  <2a> 

where  Ap  5  maximum  capillary  head 

c  ,max 

Ap  gravitational  head  acting  on  liquid  (evaporator  assumed 
9  to  be  positioned  below  the  condenser) 

Ap ^  pressure  drop  in  liquid  due  to  viscous  drag 

Ap  pressure  drop  in  vapour  due  to  viscous  drag  and  mass 

transfer . 

This  equation  states  that  for  any  fluid  path  defined  by  two 
points  along  the  heat  pipe,  the  sum  of  the  driving  forces 
(capillary  head  and  gravitational  head)  must  be  equal 
to  the  sum  of  the  pressure  losses  undergone  by  the  liquid  and 
vapour.  These  pressure  losses  are  dependent  on  the  fluid  cir¬ 
culation  rate  (and  hence  the  thermal  load),  whereas  the  maximum 
capillary  head  and  the  gravitational  head  are  essentially  in¬ 
dependent  of  load.  When  the  load  increases  so  that  the  two 
sides  of  the  equation  are  equal  anywhere  along  the  heat  pipe, 
the  capillary  pumping  limited  has  been  reached. 

Finally,  there  exists  a  boiling  (burnout)  limit  which 
corresponds  to  the  formation  of  a  stable  vapour  film  within  the 
wick  due  to  the  excessively  rapid  production  of  bubbles,  whose 
escape  is  inhibited  by  the  wick  structure.  This  condition 
tends  to  arise  for  high  operating  temperatures  where  there  is 
a  large  temperature  differential  across  the  evaporator  wick. 

In  all  cases,  the  heat  pipe  operates  within  the  area 
bounded  by  the  curves  shown  in  figure  2.  For  the  standard  heat 
pipe,  the  actual  operating  point  is  determined  solely  by  the 
heat  source  and  sink  temperatures.  The  goal  of  variable  con¬ 
ductance  designs  is  to  be  able  to  displace  this  point  at  will, 
and  in  some  cases  even  to  change  the  shape  cf  the  operating 
envelope  defined  by  these  limiting  conditions. 
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2 . 3  Variable  conductance  heat  pipes 

For  reasons  already  stated,  the  concept  of  variable 
conductance  heat  pipes  was  extensively  developed  beginning  in 
the  late  1960s.  Most  variable  conductance  designs  employ  a 
second  fluid,  an  inert  non-condensable  gas,  to  act  as  a  buffer 
to  the  heat  transfer  process.  The  inert  gas,  which  in  prin¬ 
ciple  remains  unmixed  with  the  working  fluid,  is  compressed 
into  a  reservoir  at  the  condenser  end  of  the  heat  pipe  by  the 
total  pressure  of  the  working  vapour  (Fig.  3).  The  non¬ 
condensable  gas  blocks  a  portion  of  the  condenser  whose  length 
depends  on  the  thermal  load  on  the  heat  pipe.  As  the  thermal 
load  increases,  the  inert  gas  is  further  compressed,  exposing 
more  of  the  condenser  to  the  working  vapour,  thus  increasing 
the  heat  transfer  rate  and  limiting  variations  in  heat  source 
and/or  sink  temperature. 

The  heat  transfer  rate  may  either  be  controlled 
passively,  as  just  described,  or  actively,  in  which  case  the 
inert  gas  reservoir  is  heated  or  cooled  in  order  to  regulate 
more  precisely  the  position  of  the  vapour/gas  interface  and 
hence  the  effective  length  of  the  condenser  section. 

There  exist  other  schemes  for  variable  conductance 
heat  pipes;  these  are  mostly  the  so  called  vapour  flow  modula¬ 
tion  or  liquid  flow  modulation  types,  where  a  throttling 
mechanism  is  placed  either  in  the  vapour  or  liquid  phase  of 
the  flow  circuit.  This  may  consist  of  a  variable  orifice 
through  which  the  vapour  passes  on  its  way  from  the  evaporator 
to  the  condenser,  or  it  may  be  some  sort  of  controllable 
interruption  or  impediment  to  condensate  flow  in  the  wick. 
Unlike  the  non-condensable  gas  reservoir  type  of  VCHP,  both 
of  these  latter  concepts  require  the  use  of  moving  parts. 

All  of  the  above  mentioned  methods  of  heat  pipe 
control  have  several  properties  in  common.  First,  the 
quantity  of  working  fluid  present  in  the  heat  pipe  remains 


constant  regardless  of  the  control  setting.  Secondly,  these 
control  systems  can  only  attenuate  the  heat  flux  from  some 
maximum  value  attainable  were  the  control  system  not  present. 
These  systems,  then,  meet  the  needs  of  only  one  class  of 
applications  discussed  in  the  introduction,  namely  temperature 
control,  and  this  only  when  the  thermal  potential  placed 
across  the  heat  pipes  is  always  greater  than  or  equal  to  the 
required  heat  transfer. 


2 . 4  Liquid  reservoir  variable 
conductance  heat  pipes 

The  liquid  reservoir  VCHP  is  based  on  a  different 
principle  than  the  proceeding  types.  It  uses  no  direct  inter¬ 
vention  in  the  fluid  circulation  process  to  control  heat  flux; 
instead,  the  thermal  conductance  is  varied  by  modulating  the 
heat  pipe's  vapour  temperature.  Refering  to  the  cutaway  view 
in  figure  4,  one  sees  that  the  evaporator  end  of  each  heat 
pipe  is  connected  to  an  unwicked  reservoir  capable  of  holding 
condensed  working  fluid.  Vapour  can  flow  freely  from  the 
reservoir  to  the  heat  pipe,  or  vice  versa,  whereas  liquid 
condensate  can  only  run  down  from  the  heat  pipe  into  the 
reservoir,  since  the  heat  pipe  is  operated  in  the  vertical 
position.  The  heat  pipe/reservoir  combination  contains  slightly 
more  working  fluid  than  is  required  to  completely  saturate  the 
wick;  thus,  a  liquid/vapour  interface  always  exists  within  the 
reservoir.  A  possible  exception  is  when  the  heat  pipe  is 
functioning  at  a  high  heat  transfer  rate  with  a  high  reservoir 
temperature,  a  situation  in  which  reservoir  dryout  could  occur. 
Since  heat  pipe  operation  depends  on  the  existence  of  pressure 
differentials  between  the  liquid  and  vapour  phases  in  the  eva¬ 
porator  and  condensor,  as  well  as  pressure  gradients  within 
each  phase  (as  expressed  by  equation  2.1),  it  cannot  truly  be 
said  that  the  liquid  and  vapour  phases  are  in  equilibrium. 

The  variations  in  vapour  pressure  and  temperature  are  so  small, 
however,  that  the  vapour  can  be  considered  saturated  everywhere 
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in  the  heat  pipe  and  in  the  reservoir.  Consequently,  the 
operating  temperature  (and  hence  the  operating  pressure)  is 
determined  by  the  temperature  of  the  liquid  in  the  reservoir. 
For  the  experimental  work  described  in  this  report,  the  assump¬ 
tion  was  made  that  due  to  the  high  conductivity  of  the  copper 
reservoir  wal 1 s ,  the  externally  measured  reservoir  temperature 
is  approximately  equal  to  the  vapour  temperature  in  the  heat 
pipe. 


Variable  conductance  is  achieved  in  this  context  by 
variation  of  the  reservoir  temperature,  which,  by  imposing  an 
operating  temperature  and  pressure  on  the  heat  pipe,  causes  a 
particular  operating  regime  to  be  established  by  the  interplay 
of  thermal  and  hydrodynamic  processes  in  the  heat  pipe.  This 
operating  regime,  with  its  corresponding  thermal  conductivity, 
is  characterized  by  the  amount  of  fluid  present  in  the  heat 
pipe  and  by  its  circulation  rate. 

2 . 5  The  prototype 

The  heat  exchanger  being  investigated  is  shown 
schematically  in  figure  3  and  its  photograph  appears  in 
figure  4.  It  consists  of  an  inline  array  of  sixteen  finned 
heat  pipes,  where  each  row  of  four  heat  pipes,  transverse  to  the 
flow  direction,  has  a  single  reservoir  serving  all  of  the  heat 
pipes  in  that  row.  There  is  a  divider  plate  located  midway 
along  the  heat  pipe  which  divides  the  heat  exchanger  into 
evaporator  and  condenser  sections.  As  seen  in  the  photograph, 
a  short  length  of  the  divider  plate  is  bent  over  at  the  ends, 
forming  tabs  which  effectively  block  part  of  the  heat  exchanger 
face,  thereby  creating  a  short  adiabatic  zone  between  the 
evaporator  (bottom)  and  condenser  (top).  A  summary  of  the  key 
heat  exchanger  dimensions,  as  well  as  the  known  and  estimated 
heat  pipe  parameters  is  given  in  Appendices  A  and  C. 


3.  THE  EXPERIMENTAL  INVESTIGATION 


3 . 1  The  test  faci 1 i ty 

To  make  optimum  use  of  the  VKI  laboratory  facilities, 
the  heat  pipe  heat  exchanger  was  adapted  to  the  Institute's 
already  existing  solar  energy  test  set-up  (Ref.  1).  The  modi¬ 
fied  system,  with  the  heat  exchanger  in  place,  is  sketched  in 
figure  6.  The  system  consists  simply  of  two  open  circuits, 
the  hot  air  path  passing  through  the  evaporator  side  of  the 
heat  exchanger,  the  cold  air  path  through  the  condenser  side. 
Each  circuit  is  equipped  with  a  fixed  speed  blower,  an  orifice 
plate  flow  meter,  and  a  valve  to  regulate  the  mass  flow  rate. 
The  hot  air  circuit  is  equipped  with  an  electrical  resistance 
heater  that  has  power  settings  of  2,  4  and  6  Kw.  The  cold 
air  circuit  is  equipped  with  a  water-to-air  heat  exchanger 
capable  of  cooling  the  air  to  approximately  15°C  (depending 
on  the  water  flow  rate)  before  it  passes  through  the  condenser. 
All  of  the  tests  referred  to  in  this  report  were  conducted  in 
the  counterflow  configuration. 

As  indicated  in  figure  6,  the  test  facility  is  out¬ 
fitted  with  a  microprocessor-controlled  data  acquisition 
system  linked  to  a  mi ni -computer  capable  of  doing  real  time 
data  reduction.  Both  the  microprocessor  and  the  mi ni -computer 
control  programs  (Ref.  2)  are  adaptations  of  the  programs  used 
to  perform  data  acquisition  in  the  SUN  experiments.  The 
versions  used  for  this  study  permit  the  scanning  of  up  to  48 
thermocouples  as  well  as  two  pressure  transducers,  with  a 
total  acquisition  time  of  about  30  seconds  and  a  variable 
interval  between  acquisitions. 

Temperature  control  of  the  liquid  reservoirs  was 
achieved  in  one  of  three  different  ways,  depending  on  the 
type  of  test  : 
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1.  All  four  reservoirs  were  placed  in  a  large  oil  bath 
(Fig.  7),  whose  temperature  was  maintained  nearly  uniform  by 
means  of  a  rheostaticall y  controlled  electrical  resistance 
heater  and  a  rotating  agitator. 

2.  All  four  reservoirs  were  placed  in  insulating  sleeves  to 
effectively  eliminate  heat  transfer  between  the  reservoirs 
and  the  surroundings. 

3.  The  reservoirs  were  placed  in  individual  baths,  one  of 
which  contained  heated  oil  as  in  case  1,  while  the  other 
three  contained  ice  water.  This  permitting  the  testing  of 
the  heat  exchanger  with  only  one  row  of  heat  pipes  active 
at  a  time. 


3 . 2  Instrumentation  of  the  heat  exchanger 


The  heat  exchanger  was  sandwiched  between  two  coup¬ 
ling  sections,  which  were  designed  to  provide  a  smooth  transi¬ 
tion  of  the  flow  in  the  160  mm  diameter  circular  ducting  to 
the  rectangular  flow  section  of  the  heat  exchanger.  Near  the 
inlet  and  outlet  of  each  side  of  the  heat  exchanger  there  were 
apertures  in  the  coupling  sections  which  gave  access  to  the 
test  sections.  The  whole  assembly  consisting  of  heat  exchanger 
and  coupling  sections  was  mounted  on  a  support  designed  for 
testing  the  heat  exchanger  at  0°,  30°, 45°,  60°  and  90°  from 
the  vertical.  All  measurements  discussed  in  this  report  were 
performed  with  the  heat  exchanger  in  the  vertical  position. 

Static  pressure  taps  were  drilled  in  the  walls  of 
the  coupling  sections  opposite  the  test  section  apertures 
(Fig.  7).  A  pressure  transducer  was  connected  to  each  pair 
of  taps  to  measure  the  pressure  differential  across  the  eva¬ 
porator  and  condenser  sections  of  the  heat  exchanger.  Although 

several  types  (Setra,  Validyne)  of  transducers  were  used  during 
the  courseof  experiments,  they  all  had  a  high  sensitivity 
(range  ^  70  mm  H20)  and  they  all  gavean  output  on  the  order 
of  one  volt,  so  as  to  be  compatible  with  the  data  acquisition 
system  without  the  use  of  intermediate  amplifiers. 


sm 
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Two  types  of  probes  were  placed  in  the  test  sections, 
a  vane  anemometer  to  measure  the  velocity  of  the  air  stream 
and  thermocouples  to  measure  its  temperature.  A  plexiglas 
sliding  mechanism  (Fig.  B)  was  designed  to  hold  the  anemometer 
and  seal  the  test  section  apertures,  while  allowing  manual 
displacement  of  the  probe  in  the  y  and  z  directions  (x  being 
the  flow  direction).  All  velocity  measurements  were  made 
without  heating  or  cooling  or  the  airstream  and  the  data 
acquisition  system  was  not  used.  For  the  thermal  measurements, 
thermocouples  were  placed  at  the  inlet  and  outlet  of  both 
sides  of  the  heat  exchanger.  A  plate,  through  which  were 
passed  two  thermocouple  wires,  was  placed  over  the  aperture 
at  each  inlet  section.  These  thermocouple  wires  were  drawn 
taut  by  fine  copper  wires  which  passed  through  the  static 
pressure  taps  and  were  clamped  there  by  the  plastic  pressure 
tubing.  The  wires  were  positioned  such  that  the  temperature 
sensitive  thermocouple  junctions  were  situated  midway  across 
the  test  section.  This  simple  but  effective  system  was  adap¬ 
ted  in  order  to  minimize  perturbation  of  the  flow  upstream 
of  the  heat  exchanger. 

In  both  the  evaporator  and  condenser  outlet  test 
sections  was  placed  a  grid  of  fifteen  thermocouples,  arranged 
in  five  rows  of  three  thermocouples  (Fig.  9).  The  thermocouple 
heads  protrude  from  hollow  steel  tubing,  which  can  slide  with 
respect  to  the  mounting  plate,  thus  permitting  a  degree  of 
freedom  in  the  lateral  position  of  the  thermocouple  columns. 

For  the  tests  described  in  this  paper,  the  three  columns  (in 
each  grid)  of  thermocouples  were  aligned  with  the  spaces 
between  the  columns  of  heat  pipes. 

There  were  additional  thermocouples  connected  to 
the  data  acquisition  system.  One  of  these  was,  of  course, 
a  reference,  which  was  placed  in  an  ice  bath;  the  others 
measured  the  ambient  air  temperature  and  the  external  tem¬ 
perature  of  each  reservoir. 
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3 . 3  Preliminary  measurements  : 
flow  characteri sti cs 

3.3.1  Xel°ci  ty_p£of  iHe 

Velocity  profile  measurements  were  conducted  for  two 
reasons  :  first,  to  get  an  idea  of  the  uniformity  of  the  flow 
at  the  inlet  and  outlet  of  the  heat  exchangers,  and  second, 
to  determine  the  optimal  placement  of  the  thermocouples  at 
the  outlet  section. 

It  was  found  that  when  the  probe  was  placed  directly 
downstream  of  a  column  of  heat  pipes,  the  measured  air  velo¬ 
city  varied  quite  erratically  as  the  probe  was  displaced  in 
the  z  direction,  indeed,  falling  to  zero  (or  at  least  below 
the  sensitivity  threshold  of  the  anemometer)  at  some  points 
along  the  measurement  line.  We  concluded  from  this  that 
there  remains  a  turbulent  wake  downstream  of  the  heat  exchan¬ 
ger  behind  each  column  of  heat  pipes  in  which  the  vane  ane¬ 
mometer  is  unsuitable  for  velocity  measurement.  The  non¬ 
uniformity  of  this  wake  (as  indicated  by  the  vane  anemometer) 
may  be  due  to  a  high  degree  of  sensi ti vi ty,  in  the  wake  region, 
to  small  irregularities  in  fin  geometry  and  spacing.  This 
possibility  is  suggested  by  the  fact  that  the  lightweight 
fins  had  been  deformed  at  various  points  on  both  faces  of 
the  heat  exchanger  due  to  previous  mishandling.  Due  to  the 
presence  of  these  highly  irregular  wake  regions,  it  was 
decided  to  represent  the  velocity  profiles  in  the  outlet 
section  in  terms  of  measurements  made  along  lines  in  the 
streamplanes  midway  between  the  heat  pipe  columns.  This  was 
also  the  reason  for  positioning  the  thermocouple  columns  mid¬ 
way  between  the  heat  pipe  wake  regions. 

The  two  dimensional  inlet  and  outlet  profiles  were 
measured  for  the  evaporator  section  with  a  mass  flow  rate  of 
approximately  .1  kg/s.  At  the  inlet  section  there  were 
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velocity  gradients  in  both  the  vertical  and  lateral  directions. 
A  careful  examination  of  the  heat  exchanger  assembly  revealed 
that  the  non-uniformity  of  the  inlet  velocity  profile  corres¬ 
ponded  to  a  misalignment  of  the  air  inlet  pipe.  The  outlet 
velocity  profile,  though  based  on  a  much  sparser  data  field, 
indicates  that  as  the  air  passes  through  the  heat  exchanger, 
a  high  degree  of  flow  uni formi zation  occurs  in  the  lateral 
direction,  whereas  the  average  vertical  velocity  gradients 
tended  to  remain  constant.  This  is  intuitively  comprehen¬ 
sible;  air  that  enters  fins  with  velocity  gradients  in  both 
the  y  and  z  directions  has  a  chance  to  redistribute  itself 
in  a  plane  parallel  to  the  fins,  thereby  reducing  lateral 
gradients,  but  not  in  a  plane  normal  to  the  fins,  once  the 
air  is  trapped  between  the  fins.  Hence,  except  for  some 
slight  vertical  flow  redistribution  that  may  take  place  im¬ 
mediately  upstream  of  the  heat  exchanger,  the  vertical  velo¬ 
city  gradients  are  preserved.  These  findings  lead  one  to 
conclude  that,  although  it  may  be  difficult  to  achieve  a 
totally  uniform  flow  field  at  the  inlet,  special  attention 
should  be  paid  toward  minimizing  velocity  gradients  normal 
to  the  fins. 

The  most  important  conclusion  to  be  drawn  from  these 
tests  is  that  the  outlet  velocity  profile  was  sufficiently 
uniform  (out  of  the  wake  zone)  to  justify  a  simple  averaging 
of  all  the  temperatures  given  by  the  thermocouple  grid  to 
determine  the  bulk  temperature  of  the  outlet  air  stream. 

3.3.2  P.rf  ssure_drop 

As  part  of  the  preliminary  testing  program,  a  series 
of  measurements  was  conducted  to  determine  the  functional 
relationship  between  the  mass  flow  rate  through  the  heat  ex¬ 
changer  and  the  resulting  pressure  drop.  This  information  is 
used  for  two  purposes  : 
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1.  During  the  thermal  measurements,  the  mass  flow  rates  are 
computed  from  the  pressure  drop  across  each  side  of  the  heat 
exchanger.  This  is  a  very  direct  method  of  flow  metering, 
one  which  eliminates  errors  due  to  leaks  at  other  points  in 
the  flow  circuit,  which  in  fact  were  found  to  be  a  problem 
in  the  SUN  installation. 

2.  A  derivation  of  the  convective  heat  transfer  coefficient 
in  the  heat  exchanger,  necessary  for  the  computer  model,  is 
possible  knowing  the  function  relating  mass  flow  rate  to 
pressure  drop. 

We  define  the  friction  factor  with  the  following 
equation,  which  relates  the  average  component  in  the  flow 
direction  of  the  wall  shear  stress  in  the  heat  exchanger  to 
the  dynamic  pressure  of  the  air  based  on  its  bulk  velocity 
in  the  inlet  section 

t  =  f  (3.1) 

2 

where  f  =  friction  factor. 

t  =  average  axial  component  of  wall  shear  stress  in  the 
heat  exchanger 

p  =  ai r  density 

u  =  bulk  velocity  of  air  in  inlet  section 

The  mass  flow  rate  through  the  heat  exchanger  is 

given  by  : 

m  =  PuAn  (3.2) 

where  An  =  cross  sectional  flow  area  (inlet  section) 


Taking  the  heat  exchanger  itself  as  a  control  volume, 
a  force  balance  on  the  air  gives 

4pAn  =  t  Ahx  (3.3) 

where  A^x  =  heat  exchange  surface  area. 

Combining  equations  3.1  and  3.3,  one  gets 


where  c  =  — .  Substituting  this  into  equation  3.2  gives 
AHX 


The  lengths  of  the  evaporator  and  condenser  sections  are  not 
exactly  equal;  there  are  different  values  for  Ap  and  A^x 


A  =  height  x  width 
n 

(Anj  =  .24  .13  =  .0312  m 
[Anj  =  .22  .13  =  .8286  m 

Ahx  =  2  (number  of  fins)  (exposed  fin  area)  + 

+  (number  of  heat  pipes)  (exposed  pipe  area) 
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i  =  height  of  section 

s  =  fin  spacing 

t  =  fin  thickness 

Af  =  projected  area  of  fin 

Np  =  number  of  heat  pipes 

d  =  diameter  of  heat  pipes 
P 

(ahx ) e  =  4,23  m2 

(ahx)c  -  3-88  "2 

A 

Hence,  C  =  -  =  .00738  for  both  sides  of  the  heat  exchanger. 

AHX 

Figures  11  and  12  show  the  mass  flow  rate  (as 
measured  with  a  reference  orifice  plate  flow  meter,  Ref.  3) 
plotted  with  respect  to  the  square  root  of  the  pressure  drop 
measured  across  the  evaporator  and  condenser  sections  of  the 
heat  exchanger,  respectively.  These  measurements  showed 
very  little  random  error  and  were  highly  repeatable.  Clearly, 
the  most  striking  feature  of  the  data  is  its  linearity  above 
a  certain  threshold  mass  flow  rate  of  about  .04  kg/s.  This 
means  that  in  equation  3.4,  the  friction  factor  f  can  be 
considered  constant  in  this  range,  provided  a  constant  term 
is  added  to  the  right  hand  side  of  the  equation  to  compensate 
for  the  fact  that  the  extrapolation  of  the  linear  portion  of 
the  curve  does  not  pass  through  the  origin.  There  is  no 
a  p> iori  reason  that  this  must  be  so;  indeed,  the  friction 
factor  in  any  internal  flow  situation  is  a  function  of  Reynolds 
number.  It  appears,  however,  that  above  a  certain  mass  flow 
rate,  the  friction  factor  for  this  particular  geometry  achieves 
Reynolds  number  independence,  as  does  the  friction  factor  for 
turbulent  flow  in  rough  pipes. 


A  best  fit  to  the  points  in  the  linedr  range  gives 
the  following  expressions  for  mass  flow  rate. 

Condenser 


m  =  .0142  \/Tp  -  .0036  [  a  pi  =  Pa 

™  =  .0444  vTp  -  .0336  [  a p]  =  mm  H2Q 

m  =  .  2 0 2 £  c  >/^p  -  .0336  [  a p]  =  mm  H20 


(3.5a) 


(3.5b) 


(3.5c) 


[  £ C1  =  meters 


Evaporator 


m  =  .0150  v/Tp  -  .0340  [  a p]  =  Pa 

m  =  .0471  vTP  -  .0340  [  a p]  *  mm  H20 

m  =  .  0 1 96 i c  VIp  -  .0340  [  a p]  =  mm  H20 


(3.6a) 


(3.6b) 


(3.5c) 


[  £e]  *  meters 

The  slope  of  these  lines  is  equal  to  the  grouping 
Anv/20c/f ,  from  which  f  can  be  determined.  The  following 
values  are  found. 

condenser  :  f  =  .078 
evaporator  :  f  =  .074 


Since  the  friction  factor  (assuming  Reynolds 
independence)  is  only  a  function  of  flow  geometry,  it  is  rea¬ 
sonable  to  assume  that  the  friction  factor  for  both  sides  of 
the  heat  exchanger  should  be  the  same.  For  simplicity,  then, 
we  will  use  the  average  value  f  =  .076. 


3.4  Thermal  measurements 


3.4.1  Da  ta_a£qu_i  si;t  i£n_a  nd_proce£si_n£ 

Figure  12  is  a  close-up  view  of  the  data  acquisition 
system  (DAS),  which  rapidly  gathers  and  processes  all  relevant 
pressure  and  temperature  data.  Figure  13  shows  the  heat 
exchanger  mounted  on  its  support,  outfitted  with  thermocouples 
and  pressure  transducers,  and  connected  to  the  DAS.  During 
each  acquisition  cycle,  the  ME  K 1 1  microprocessor  feeds  data  to 
the  PET  minicomputer,  where  BASIC  control  program  then  performs 
calculations  and  prints  out  the  results.  What  follows  is  a 
description  of  the  calculations  in  the  data  reduction  process. 

Mass  flow  rate 


The  mass  flow  rates  through  the  evaporator  and 
condenser  sections  are  computed  on  the  basis  of  the  pressure 
transducer  output  using  equations  3.5c  and  3.6c. 

Heat  transfer  rate 


The  heat  transfer  rate  between  air  streams  is  not 
measured  directly.  The  measurements  only  allow  the  calcula¬ 
tion  of  the  heat  given  up  by  the  hot  air  and  the  heat  absorbed 
by  the  cold  air.  Since  the  heat  pipes  are  connected  to  liquid 
reservoirs,  which  may  act  as  either  heat  sources  or  sinks, 
these  two  heat  transfer  rates  are  not  necessarily  the  same. 
They  are  given  by  the  expressions 


Qe 

=  m  c  | 
e  p  1 

(T:,e-T2,e) 

• 

=  me  1 
c  p  1 

iT2.c*Tl.c) 

(3.7) 

(3.8) 


The  bulk  outlet  temperatures  T,  and  T,  _  are  determined  as 

t  )  c  t  j  C 

mentioned  in  section  3.3.1,  by  averaging  the  15  values  given 
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by  each  thermocouple  grid.  The  reservoir  power  is  defined. as 
the  difference  between  the  condenser  and  evaporator  heat 
transfer  rates  : 

Qr  -  Qc-Qe  (3.9) 


Effectiveness 


The  effectiveness  of  a  heat  exchanger  is  a  dimension¬ 
less  parameter  used  to  compare  the  performance  of  different 
heat  exchangers  or  of  the  same  heat  exchanger  under  different 
operating  conditions.  This  quantity,  denoted  e,  is  defined 
as  the  ratio  of  the  actual  heat  transfer  rate  in  the  heat 
exchanger  to  the  maximum  thermodynamically  possible  heat 
transfer  rate.  The  maximum  possible  heat  transfer  rate  is  given 
by  the  maximum  temperature  difference  existing  within  the  heat 
exchanger,  the  difference  between  the  two  inlet  temperatures, 
times  the  lower  value  of  the  heat  capacity  rate  of  the  flow  : 


In  the  present  case,  the  heat  capacity  Cp  of  the  air  is 
assumed  constant;  thus,  (mcp)m^  n  always  pertains  to  the  flow 
with  the  minimum  mass  flow  rate. 


The  two  principal  heat  transfer  rates  in  the  heat 
exchanger,  Q  and  Q  ,  lead  to  two  different  definitions  of 

C 

the  effectiveness  : 

Qe 

Evaporator  effectiveness  e  -  - 


Condenser  effectiveness 


'max 


There  are  two  dimensionless  parameters  computed 
during  the  data  reduction  process  which  are  useful  both  in 
prescribing  the  test  conditions  and  organizing  the  effectiveness 
data.  There  are  the  mass  flow  ratio  and  the  normalized  reser¬ 
voir  temperature. 

Mass  flow  ratio,  MR 

The  mass  flow  ratio  is  defined  as  the  air  mass  flow 
rate  through  the  condenser  section  divided  by  the  air  mass  flow 
rate  through  the  evaporator  : 

m 

MR  = 


Normalized  reservoir  temperature,  Tn 


Since  the  reltionship  of  reservoir  temperature  to 
effectiveness  is  of  primary  interest  in  this  study,  it  is 
necessary  to  express  reservoir  temperature  in  some  way  that 
permits  comparisons  of  tests  under  different  thermal  loads  for 
various  reservoir  temperatu res .  In  order  that  the  chosen 
transformation  be  valid,  it  must  be  independent  of  the  perfor¬ 
mance  of  the  heat  exchanger,  i.e.,  it  must  be  independent  of 
all  temperatures  downstream  of  the  inlets.  This  leaves  little 
choice  but  to  define  a  parameter  that  indicates  the  reservoir 
temperature's  relative  position  between  the  cold  and  hot  air 
inlet  temperatures  : 


VTl,c 

T1 ,e"Tl ,c 


It  should  be  mentioned  here  that  the  liquid  reservoir 
VCHP  heat  exchanger  could  be  considered  a  reservoir  power 
modulated  system  rather  than  a  reservoir  temperature  modulated 


one,  in  which  case  it  would  be  appropriate  to  define  a  norma¬ 
lized  reservoir  power,  Qp/Qmax*  A  somewhat  different  test 
facility  would  have  been  required  to  conduct  measurements  from 
this  point  of  view,  one  in  which  reservoir  power  could  be 
accurately  measured  (now  it  can  only  be  determined  by  taking 
the  difference  between  Qc  and  Qg,  which  already  are  of  only 
marginal  accuracy).  The  author  feels  that  the  temperature 
formulation  is  more  fundamental  to  the  physical  processes 
responsible  for  the  variable  conductance. 

3.4.2  GJoba^Jiea t  jexcha_n£er  £erforma_nce  - 
Xestjs  i_n_l ar£e_o_U_b£th^ 

A  test  matrix  was  established  in  which  each  point 
was  defined  by  three  main  parameters  :  heater  power,  mass  flow 
ratio  and  reservoir  temperature.  The  heat  exchanger  was  tested 
at  thermal  loads  correspondi ng  to  heater  settings  of  2,  4  and 
6  kW  :  at  mass  flow  ratios  of  .4,  .6,  .8,  and  1.0;  and  at 
irregular  intervals  in  the  normalized  temperature  range 
0  <  Tn  <  1 .  In  addition,  several  tests  were  run  at  Tn  <  0  and 
Tp  >  1  to  see  the  effect  of  operating  outside  the  normal  reser¬ 
voir  temperature  range. 

For  each  test  point,  the  data  acquisition  system  was 
activated  once  the  test  conditions  (valve,  air  heater  and 
rheostat  settings)  were  established.  The  program  output  (heat 
transfer  rates,  effectiveness  and  temperature  profile)  was 
monitored  until  steady  state  was  reached,  typically  taking 
about  20  minutes;  then  the  DAS  was  stopped  and  reactivated, 
with  subsequent  acquisitions  being  recorded  on  magnetic  cas¬ 
sette.  An  example  af  the  DAS  output  appears  in  Appendix  B. 

The  evaporator  and  condenser  effectiveness  for  MR  =  1 
are  plotted  in  figure  14.  The  points,  although  correspondi ng 
to  tests  under  different  thermal  loads,  admit  of  a  relatively 
smooth  curve  fit,  which  validates  the  method  of  non-dimensiona- 
lizating  the  reservoir  temperature.  The  two  curves,  evaporator 
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and  condser  effectiveness,  cross  at  a  point  designated  the  • 
adiabatic  operating  point.  At  normalized  reservoir  tempera¬ 
tures  below  the  adiabatic  operating  temperature,  the  evaporater 
effectiveness  is  higher  than  the  condenser  effectiveness,  the 
difference  corresponding  to  the  rate  of  heat  transfer  from  the 
heat  pipes  to  the  reservoirs.  Under  steady  state  conditions, 
this  heat  loss  can  only  be  by  conduction,  since  if  it  were  by 
convection  the  heat  pipe  wick  would  gradually  empty  of  liquid. 
This  corresponds  to  a  condition  in  which  there  is  dryout  at  the 
lower  end  of  the  evaporator  section  in  at  least  some  of  the  heat 
pipes,  a  condition  manifested  by  a  steep  gradient  in  the  air  tem¬ 
perature  profile  measured  at  the  exit  of  the  evaporator  section 
(Fig.  15).  The  temperature  profile  shown  probably  represents 
a  superposition  of  several  staggered  temperature  steps,  since 
the  dryout  boundary  (if  any)  may  be  at  a  different  location  in 
each  row  of  heat  pipes. 

At  the  adiabatic  operating  point,  there  must  be  no 
influx  of  liquid  into  the  reservoir,  since  that  would  neces¬ 
sitate  reservoir  heating  in  order  to  revaporize  the  fallen 
liquid.  There  are,  then,  two  possible  adiabatic  operating 
situations  :  One  is  where  the  reservoirs  are  completely  dried 
out,  with  all  of  the  liquid  circulating  in  the  heat  pipes. 

In  this  case,  heat  input  to  the  reservoir  serves  only  to  main¬ 
tain  it  at  an  elevated  temperature  (near  the  hot  air  inlet 
temperature),  compensating  for  heat  losses  by  conduction 
radiation,  etc.  The  second  is  where  there  is  some  condenser 
liquid  in  the  reservoir;  the  adiabatic  operating  temperature 
in  this  case  is  that  at  which  the  reservoir  liquid  is  in 
equilibrium  with  the  heat  pipe  vapour.  This  situation  is 
plausible  with  or  without  dryout  in  the  heat  pipe  evaporator 
section . 


Above  the  adiabatic  operating  point,  the  condenser 
heat  transfer  is  greater  than  that  of  the  evaporator,  being  the 
sum  of  the  heat  input  to  the  evaporator  and  the  heat  input  to 
the  reservoirs.  The  true  effectiveness  of  the  heat  exchanger, 
at  any  given  operating  temperature,  must  be  regarded  as  the 
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lower  of  these  two  curves,  since  it  is  always  the  lower  curve 
that  represents  the  heat  acutally  exchanged  between  the  two 
air  streams,  not  taking  into  account  heat  which  might  have 
been  lost  or  gained  from  the  reservoirs.  The  significant 
feature  of  the  super-adiabatic  operating  range  (T  >  Ta)  is 
that  after  increasing  nearly  linearly  with  T  ,  the  evaporator 
effectiveness  reaches  a  peak  and  then  decreases  as  the  norma¬ 
lized  reservoir  temperature  approaches  1.0.  This  decline  in 
evaporator  effectiveness  is  explained  by  the  fact  that  at 
high  normalized  reservoir  temperatures  there  is  a  low  thermal 
gradient  across  the  liquid  layer  in  the  evaporator  section, 
and  hence  a  low  heat  flux.  The  maximum  in  the  evaporator 
effectiveness  curve  indicates  the  maximum  heat  exchange  rate 
between  air  streams  that  can  be  achieved  by  the  heat  exchanger 
for  these  inlet  conditions  (and  with  uniform  reservoir 
temperature) . 

The  usable  heat  flux  range  can  effectively  be  exten¬ 
ded  if  one's  principal  purpose  is  to  heat  cold  air.  Although 
the  maximum  heat  recovery  from  the  hot  air  is  given  by  the 
evaporator  effectiveness,  the  heat  supplied  to  the  reservoirs 
is  also  recovered  at  the  condensor;  thus,  in  cases  where  the 
thermal  potential  of  the  evaporator  is  insufficient  to  meet 
the  heat  demand  at  the  condenser,  power  can  be  supplied  to 
meet  the  deficit.  It  should  be  noted,  however,  that  this  is 
an  inefficient  mode  of  operation,  since  increasing  the  norma¬ 
lized  reservoir  temperature  beyond  Tmax  reduces  the  effective¬ 
ness  of  the  evaporator,  so  the  additional  reservoir  power 
required  to  attain  the  desired  rate  of  heat  transfer  at  the 
condenser  will  be  greater  than  the  original  deficit. 

Figure  16  is  a  plot  of  the  true  heat  exchanger 
effectiveness  versus  normalized  reservoir  temperature.  It 
represents  the  ensemble  of  the  large  oil  bath  test  results; 
all  thermal  loads  and  mass  flow  ratios  are  included.  The 
effect  of  decreasing  the  ratio  of  mass  flow  rates  from  1.0 
is  to  displace  the  effectiveness  curve  upward  and  slightly 
to  the  right,  giving  both  a  higher  maximum  effectiveness  and 
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higher  optimal  operating  temperature.  The  rise  in  effectiveness 
with  decreasing  mass  flow  ratio  is  typical  of  all  fluid-to- 
fluid  heat  exchangers. 

3.4.3  Ad_i_abat;i c  op^rati ng  point  -  Individuality 
i_n£uj_aj:e^  reservoirs 

For  these  tests  each  reservoir  was  wrapped  in  a 
tightly  fitting  sleeve  made  of  closed-cell  foam.  This  indepen¬ 
dent  insulation  allowed  each  row  of  heat  pipes  to  establish  its 
own  adiabatic  operating  regime.  Indeed,  it  was  found  during 
these  tests  that  each  heat  pipe  row  had  a  slightly  different 
adiabatic  operating  temperature.  To  compare  the  results 
obtained  under  these  adiabatic  test  conditions  with  the  "adia¬ 
batic  points"  given  by  the  intersection  of  the  evaporator  and 
condenser  effectiveness  curves,  one  must  average  the  tempera¬ 
tures  of  the  insulated  reservoirs  in  order  to  obtain  an 
equivalent  normalized  reservoir  temperature.  The  comparison 
is  therefore  not  rigrous,  and  some  data  suggests  that  it  may 
be  misleading  in  the  case  where  the  heat  transfer  rate  is  a 
small  fraction  of  the  heat  exchanger's  rated  capacity.  The 
point  corresponding  to  the  highest  heat  transfer  rate  tested 
(MR  =  1,  air  heater  power  =  6  kw)  is  plotted  against  the 
general  effectiveness  curves  in  figure  14. 

An  unexpected  and  still  unexplained  finding  made 
during  the  tests  with  individually  insulated  reservoirs  is 
that  the  first  and  third  reservoirs  had  a  higher  temperature 
than  the  second  and  fourth,  the  first  row  being  the  one  closest 
to  the  hot  air  inlet.  One  would  expect  the  reservoir  tempera¬ 
tures  to  either  increase  or  decrease  monotonical ly  from  one 
row  to  the  next,  depending  on  the  direction  of  the  flow.  That 
this  is  not  the  case  indicates  that  for  some  reason  the  heat 
pipe  rows  are  not  all  operating  in  the  same  regime,  as  charac¬ 
terized  by  their  respective  normalized  reservoir  temperature 
evaluated  with  respect  to  the  local  evaporator  and  condensor 
temperatures . 
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The  results  of  these  tests  are  summarized  in  table  1. 
It  is  worth  noting  that  the  heat  losses  associated  with  adia¬ 
batic  operation  increase,  as  a  percentage  of  condenser  effec¬ 
tiveness,  as  the  overall  heat  transfer  rate  decreases.  This 
is  seen  in  the  relatively  large  spread  between  the  evaporator 
and  condenser  effectiveness  found  in  the  last  five  rows  of  the 
tabl e . 


3.4.4  Performance.  of_a_sj_n£l£  row  of_hea £  pi£es_  - 
J^n  d  i  v.i  d  u  £  1  _reservoj_r_batjis 

To  compare  the  individual  performance  of  each  row 
of  heat  pipes,  four  tall  narrow  baths  were  used  for  reservoir 
temperature  control,  three  of  them  containing  ice  water  and 
the  fourth  containing  oil  with  a  resistance  heater.  The  reser¬ 
voirs  were  tested  consecutively  under  the  following  (nominal) 
conditions  :  MR  =  .8,  heater  power  =  4  kW,  normalized  reservoir 
temperature  =  0.7. 

The  test  results  are  summarized  in  table  2.  They 
show  that  there  are  no  large-order  differences  in  effectiveness 
among  the  four  rows  of  heat  pipes.  It  is  not  clear  whether  the 
differences  that  do  exist  result  from  differences  in  heat  flux 
path  length  in  the  fins,  which  varies  with  the  position  of  the 
heat  pipe  row  tested,  or  rather  differences  in  the  heat  pipes 
themselves,  such  as  non-uniformities  in  wick  porosity,  quantity 
of  working  fluid,  or  in  other  basic  heat  pipe  characteristics. 
This  question  could  be  resolved  with  a  numerical  model  of  the 
heat  exchanger,  capable  of  simulating  single  row  operation. 
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4.  DEVELOPMENT  OF  A  NUMERICAL  MODEL 
OF  THE  HEAT  EXCHANGER 


4 . 1  Simplifying  assumptions 


In  addition  to  the  experimental  investigation 
described  above,  it  was  undertaken  to  numerically  model  the 
heat  exchanger.  This  represented  a  substantial  challenge, 
for  if  the  detailed  behaviour  of  a  single  liquid  reservoir 
VCHP  is  not  fully  understood,  the  overall  comportment  of  a 
finned  tube  heat  exchanger  containing  sixteen  such  heat  pipes 
in  cross-flow  is  understood  still  less.  Having  no  prior 
indications  as  to  the  feasibility  or  validity  of  a  complete 
computer  model  of  the  heat  exchanger,  we*  considered  it 
appropriate  to  begin  with  the  following  simplifying 
assumptions  : 

(1)  There  is  no  longitudinal  temperature  gradient  in  either 
the  evaporator  or  condenser  sections  of  the  heat  exchanger. 
Similarly,  there  are  no  transverse  temperature  gradients.  The 
only  non-uniform  temperature  profiles,  whether  of  the  fins  or 
of  the  air,  are  in  the  flow  direction. 

(2)  As  in  the  experimental  study,  only  steady  state  operation 
is  considered. 

(3)  The  liquid  and  vapour  phases  in  the  reservoirs  are  in 
equi librium. 

Consequently,  there  is  no  mass  interchange  between 
the  heat  pipes  and  the  reservoirs.  This  condition  corresponds 
to  the  adiabatic  operation  described  in  sections  3.4.2  and 
3.4.3. 


This  program  was  developed  in  conjunction  with  Paul-Herve 
Theunissen,  doctoral  candidate  at  the  Universit§  Libre  de 
Bruxelles. 
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4 . 2  Representation  of  the  fin/tube  system 

As  a  consequence  of  assumption  1  of  the  preceeding 
section  we  may  effectively  compress  all  of  the  fins  in  the 
evaporator  section  into  a  single  flat  plate,  and  likewise 
those  in  the  condenser  section.  Because  the  fins  are  thin 
and  have  high  thermal  conductivity,  the  plate  temperature 
is  assumed  to  be  constant  across  its  thickness.  These  two 
plates,  one  representing  the  evaporator  side  of  the  heat  ex¬ 
changer,  the  other  the  condenser  side,  are  coupled  to  each 
other  by  the  four  rows  of  heat  pipes.  To  facilitate  an  ana¬ 
lytical  (and  therefore  more  elegant  and  efficient)  solution 
of  the  air  and  fin  temperature  profiles,  this  coupling  is 
simulated  as  a  series  of  four  linear  heat  sinks/sources  normal 
to  the  flow  direction,  their  position  corresponding  to  that 
of  the  heat  pipe  rows.  These  four  linear  heat  sinks  end 
sources  subdivide  the  evaporator  and  condenser  plates,  res¬ 
pectively,  each  into  five  "platelets"  (Fig.  17),  for  which 
an  analytical  solution  to  the  convective  and  conductive  heat 
transfer  equations  can  be  developed. 

Because  of  the  high  thermal  conductivity  of  the  fins 
and  their  large  surface  area  relative  to  the  heat  pipe  area, 
we  assume  that  all  of  the  heat  transferred  between  the  air  and 
the  heat  pipes  passes  through  the  fins.  In  reality,  the  iso¬ 
thermal  lines  in  the  fins  (plates)  would  tend  to  form  contours 
around  the  heat  pipes;  due  to  the  linearization  of  the  heat 
si  nks/'scurces ,  however,  we  make  the  further  idealization  that 
the  isothermal  lines  are  linear  and  perpendicular  to  the  flow 
direction  (in  accordance  with  assumption  1  of  the  proceeding 
section) . 


4.2.1  The  Plate  equations 

In  any  given  platelet,  the  two  governing  heat  trans 
fer  equations  are  the  one  dimensional  conduction  equation, 
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ke  —  +  h(T  -T)  =  0 
dx2  3 

and  the  convective  heat  transfer  equation. 


me  dT 

_ £ _ i 

w  dx 


+  h(Ta-T)  =  0 


(4:1a) 


(4.2a) 


where  k  =  thermal  conductivity  of  plate 

e  =  plate  thickness  =  (thickness  of  fin)  x 
(number  of  fins) 

h  =  total  convective  heat  transfer  coefficient  = 
(convective  heat  transfer  coefficient)  x 
(number  of  fin  sides) 

w  =  width  of  plate 

m  =  total  mass  flow  rate  through  one  side  of  heat 
exchanger 

T  =  local  plate  temperature 

T  =  local  bulk  air  temperature 

a 

These  equations  can  be  simplified  using  the  following 
non-dimensional izations  : 


wh 

x  =  ~ — 
mcp 


w2  ehk 


giving 


d  2T 

“■  +  «(T  -T)  -  0 
dx2  a 


—  +  T  -T  =  0 
dx  a 


(4.1b) 


(4.2b) 
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4.2.2  So2ution_fo r_the_p2ate  temperature,  T 

i,  L. 

Consider  figure  18,  which  represents  the  itn 

platelet. 

Let  T0  «  T^O) 

Ti  =  T  (X.) 


Integrating  equation  4.1b  once,  after  substituting  the  cerm 
(Ta“T ) ,  one  obtains 


d IT 
dx 


a  ( T  -T 
v  a  ao 


) 


or 


III  + 

fT  -  I  «. 

* 

a  dX 

90  .  dx 

T 


a 


Define  T 

a 


1  dT 


then 


(4.3) 


i  d_L 

c.  dx 


+ 


T 

ex 


T 


a 
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Substituting  back  into  equation  (4.2b),  one  gets 


d  2T  +  dT 

dx2  dx 


or 


aT  +  aT  =0 

a 


d2(T-TJ  d (T-T  ) 


dx' 


dx 


a(T-T  )  =  0 
a 


Equation  4.4  is  a  second  order  linear  differential  equation 
with  the  characteristic  equation 

r2  +  r  -a  =0 


The  general  solution  of  this  equation  is 


0)1  X  o)  2  x 

T-T  =Ae  +Be 

a 


where  uj 


S  = 


to. 


■  1  +  v  1  +  4  a 


(4 . 6  a  b ) 


4.2.3  Determj nation  of_con_s tan t s_A_and  B 

We  start  with  the  implicit  boundary  conditions  on 
T0  and  Tj  : 

T0-  T  =  A+B 

a. 

1  X  oi  2  X 

T,  -  T  =  Ae  +  Be 

1  :t 


Define  3i  = 


u2X  w  i  X 

e  -e 


and  b2  = 


w2X-  w i X 

e  e 


(4. 9ab) 


then 


A}  0  i  +  A  o  $  2 


If  one  defines  A:  =  T i - T a Q  and  A0  =  T0-Ta0>  it  can  be  shov/n 
that 


dT  A i B i + A  o  &  2 

dx  ,.n  1  -  1  (ei+M 


(4.10) 


Thus , 


dx  v=0 


hence  Ta,  has  been  determined  in  terms  of  the 


"known"  temperatures  T0,  Tj  and  Tao,  thereby  enabling  a  complete 
solution  of  the  platelet  temperature  profile. 


Complete  plate  temperature  profile  given  by  equation  4.B 


34 


4.2.5  Solution  for  the  air  temperature,  T 

a 

The  expression  for  the  air  temperature  is  built  up 
from  equations  4.2b  and  4.5. 


+  T,  =  T+  Ae  +  Be 

d  Ct 


(4.11) 


The  general  solution  to  this  differential  equation 
is  of  the  form 


Ta  =  Ta  +  Ce"x  +  Cj  e  +  c2e 


(4.12) 


The  particular  solution  is  obtained  as  follows  : 


u 1 X  u2  X  Uju  w2X 

(Ci+WjCi  )e  +(^2+u)2?2)e  =  Ae  +Be 


6i  = 


e2  = 


The  constant  C  is  determined  by  setting  x  equal  to 


zero  : 


Ta0  -  T  -  C  +  ^  +  C 2 


C  =  Ta0-Ta-Ci  -  C  2 


C  .iii  -  Sl  -  S2 

a  dx 


The  air  temperature  at  the  platelet's  trailing  edge  is  then 
given  as 


+  52  e 


(4. 


Ta,  =  T  +  Ce  +  5 1  6 

1  a  1 


4.2.6  The  temperature  gradient  in  the  platelet 
at  the  trailing  edge 


Equation  4.10  expresses  the  temperature  gradient 
(which  is  proportional  to  the  heat  flux)  in  the  platelet  at 
its  leading  edge.  A  similar  expression  can  be  obtained  for 
the  trailing  edge. 


dT 


to  i  X  to 2X 

coiAe  +  to2Be 


(Oil 

U3ZX 

to2.Y 

«i*l 

(01  p 

A  i  +  A  q  e 

+  to2e 

E> 

1 

O 

o 

i _ 

to  i  X  u  2T 

Define  Sj  =  -toje  +  <o2e 


a2X  to  i  X 

e  -e 


(  a)  j.  +oj  2  )* 

(toi  -to2  )e 


to  2  X  to  i  X 

(e  -e  ) 


Then 


d I 

dx 


A i S i  +  A  q  S  2 


(4. 
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4.2.7  Ujie^nz^t^on  of_the_equations 

Equations  4.10,  4.13  and  4.14  are  all  implicitly 
linear  equations  in  terms  of  the  three  temperatures,  T0,  T,, 
and  Tao»  which  up  to  now  have  been  treated  as  known  boundary 
conditions.  The  explicit  expression  of  ^  ^  and 

in  terms  of  these  three  variable  involves  length  algebraic 
manipulations  that  will  not  be  shown  here.  The  resulting 
equations  are  given  below  : 


=  Y1T1+Y2T0"(yi+Y2)T; 


(4.15) 


where  5 


1-  -  ( B i +  S  2 ) 


Y2  = 


l-  -  ( 3 1  +  e 2 ) 


(b)  —  =  UiTi-u2To'(ui+U2)Tao 

dX  y 


(4.16) 


pi  - 


S,  + 


(Si+S2)yi 


P  2  “  S  2  + 


(Si+S2)y2 


(4-17) 


(c)  Tfli  =  v1T1+v2T0- (v1+v2-1 )Ta0 


V1  =  e3-e1+Yi£6 

v2  =  e2-e4+72e6 


el 


-x  to !  a: 

-e  +e 


0)2I 

(1+oji )  (e 


) 


0J2Z 

e2  =  e  ei 


£3 


"X  uj2X 

(~e  +e  ) _ 

io2^  w !  X 

( l  +  )  ( e  “e  ) 


£4  = 


to  i  AT 

e  e  3 


-X 

e  R  =  e  -1 


£2  +  e  3  +  e5~el “e4 

£6  =  - 

a 


4.2.8  A  s  s_emb]_y_of  the  pi  a  t£l  e 

Up  to  now,  the  entire  analytical  development  has  been 
with  regard  to  a  single  platelet,  a  segment  of  the  composite 
plate  which  represents  one  side  of  the  heat  exchanger.  These 
segments,  or  platelets,  are  joined  by  lines  of  heat  addition 
or  removal,  which  represent  the  rows  of  heat  pipes.  For  each 
side  of  the  heat  exchanger,  our  mathematical  model  links  the 
platelets  together  using  the  physically  imposed  boundary  con¬ 
ditions  of  continuity  of  temperature  in  the  air  and  conservation 


- - - 


of  energy  in  the  platelet.  At  each  node  (boundary  between, 
platelets),  the  following  conditions  apply  : 


i  +  1 


(4.18) 


kwe 


dT 


dx 


dT 

dx 


i  +  1 


X . 

1 


(4.19) 


where  =  heat  addition  rate  at  the  boundary  between  platelets 
i  and  i  +  1  (heat  transfer  rate  in  the  i^h  row  of  the 
pipes) 

The  boundary  conditions  at  the  entry  and  exit  of  the 
heat  exchanger  are  the  following: 

(1)  At  the  leading  edge,  the  plate  temperature  is  equal  to 
the  air  temperature.  This  temperature  is  not  the  upstream 
inlet  temperature  of  the  air;  due  to  stagnation  effects  at 
the  leading  edge  of  the  fins,  the  air  is  preheated  or  pre¬ 
cooled  before  actaully  passing  the  xo  =  0  line. 

(2)  At  the  trailing  edge,  it  is  assumed  that  there  is  no 
longitudinal  temperature  gradient  (no  heat  conduction)  in  the 
fin.  The  exit  air  temperature  is  given  by  an  equation  of  the 
type  4-17. 

The  above-stated  conditions  lead  to  a  system  of  twelve 
equations  in  twelve  unknowns  for  each  side  of  the  heat  exchanger. 
Referring  to  figure  17,  the  equations  for  the  top  plate  (con¬ 
denser  section)  are  as  follows  : 


Letting  $  = 


entry 


row  1 


row  2 


row  3 


kwe 


ex  i  t 


T  0 
I  n 


T° 

ao 


T!  =  T,  + 
a  o  3“ 


dT 


cx  LdX. 


<?!  - 


'dT" 

0 

rdTi 

_dX_ 

.dx. 

ao 


=  T. 


'dT' 

l 

[-dT] 

$2  ~ 

.dx. 

.dx. 

2  1 

Tao  =  Ta 


*3 


dT  1  2 

*2 


dx 


'dT' 


dx 


T  =  T 


ao  a  i 


fd  T 

3 

'it 

$4  = 

.dx.' 

l  dx. 

row  4 

*3 

t' 

a  o 


dT' 

[dx. 


‘f 

=  JviTi+v2T0-(v]+\>2_l)Tao] 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 


A  similar  set  of  equations  is  obtained  for  the  lower  plate.' 

To  solve  this  system  of  equations  a  matrix  is  composed  of  the 
coefficients  y,  y,  and  v  given  in  section  4.9.'  The  matrix 
equation  takes  the  form  shown  in  figure  19,  non-zero  elements 
indicated  by  squares  and  is  given  in  tensorial  notation  by 

A . x  =  BT  +  Ca 

—  —  —  d  0  — 

We  have,  then,  a  system  of  equations  in  which  the  only 
appearing  unknowns,  the  vector  x,  are  the  air  and  fin  tempera¬ 
tures  at  each  discontinuity  in  the  plate,  precisely  the  values 
needed  to  couple,  by  means  of  a  mathematical  model  of  the  heat 
pipes  (to  be  described  in  the  next  section),  the  two  sides  of 
the  heat  exchangers. 


4 . 3  Numerical  simulation  of  the  heat  pipes 

Ideally,  it  would  have  been  desirable  to  develop  a 
complete  numerical  model  of  the  variable  conductance  heat  pipe, 
capable  of  representing  its  operation  in  all  possible  regimes, 
thereby  permitting  a  complete  comparison  with  the  experimental 
results.  Due  to  the  complexity  of  the  internal  heat  pipe  phe¬ 
nomena,  however,  such  a  model  was  beyond  the  scope  of  this 
preliminary  study.  Instead,  we  limited  ourselves  to  a  model 
of  the  heat  pipe  in  its  adiabatic  reservoir  regime,  i.e.,  the 
point  at  which  heat  is  exchanged  only  between  the  evaporator  and 
condenser  sections.  This  approach  is  justified,  moreover,  by 
the  fact  that  we  have  been  until  now  unable  to  obtain  much  of 
the  information  on  the  physical  parameters  of  the  heat  pipe 
construction  (eg.  wick  thickness,  wick  porosity,  effective  pore 
radius,  wick  material,  etc)  necessary  to  an  accurate  prediction 
of  its  performance.  Thus,  even  restricting  ourselves  to  the 
adiabatic  operation  of  the  heat  pipes,  we  could  only  hope  to 
achieve  qualitative  agreement  with  the  experimental  results  by 
making  appropriate  guesses  at  the  unknown  parameters. 


Our  model  is  based  on  the  simultaneous  resolution 
of  the  equations  expressing  the  relevant  thermodynamic  and 
hydrodynamic  phenomena  in  the  heat  pipe.  To  attack  the  prob¬ 
lem,  it  was  first  necessary  to  have  a  geometrical  description 
of  the  liquid  layer  in  the  wick.  On  the  basis  of  a  survey  of 
the  longitudinal  temperature  distributions  we  measured  under 
adiabatic  operating  conditions,  a  number  of  different  ideali¬ 
zed  liquid  profiles  seemed  plausible  (see  Fig.  20).  Number  3 
is  the  only  one  compatible  with  the  assumption  of  constant 
heat  flux  over  the  entire  evaporator  section  and  over  the 
entire  condenser  section,  a  condition  necessary  to  the  in¬ 
corporation  of  this  model  into  the  finned  tube  model  described 
in  section  4.1.  What  follows  is  the  development  of  heat  and 
mass  transfer  equations  for  the  model  based  on  profile 

4.3.1  £rjma ry  and  s_e ccm d a_r^  n k_n o win s 

These  equations  are  formulated  in  terms  of  three 
primary  unknowns  and  three  secondary  unknowns.  The  primary 
unknowns  are 


thickness  of  the  liquid  layer  in  the  condenser  wick 
thickness  of  the  liquid  layer  in  the  evaporator  wick 
working  fluid  mass  flow  rate  through  the  adiabatic 


zone,  (m) 


max 


The  secondary  unknowns  are  the  temperatures  on  which 
the  various  fluid  properties  such  as  density  and  viscosity  are 
dependent . 


Ty  vapour  temperature 

[T  !  average  temperature  of  the  liquid  in  the  evaporator 


average  temperature  of  the  liquid  in  the  condenser 
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4.3.2  JherniaJ_  £e£i stance_network 

The  heat  transferred  by  the  heat  pipe  must  traverse 
a  network  of  thermal  resistances,  which  is  represented  in 
figure  21. 

Tg  temperature  applied  to  the  external  surface  of  the 
evaporator 

Tc  temperature  applied  to  the  external  surface  of  the 
condensor 

Tv  temperature  of  vapour  in  core  of  heat  pipe  (operating 
temperature) 

RP  thermal  resistance  across  the  heat  pipe  wall  in  the 
evaporator  section 

RW  thermal  resistance  to  conduction  across  the  liquid  layer 
in  the  evaporator  section 

RPc  thermal  resistance  across  the  heat  pipe  wall  in  the 
condensor  section 

R  thermal  resistance  to  axial  conduction  in  the  wall  of 
c 

the  heat  pipe 

For  a  heat  pipe  in  which  there  is  not  significant 
dry  out  the  axial  conduction  in  the  walls  is  negligible  .with 
respect  to  the  heat  transfer  by  vapour  transport,  so  Rc  is 
neglected.  Thus 

T  -T 

Q  =  mL  =  - - — - -  (4.32) 

RPe+RWe+RPc+RWc 

where  the  four  resistances  are  given  by 

en(Ro/Rp) 

Rp  =  - 

2*Ve 
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RP 


en( Ro/Rp) 
2lT'<p£c 


in 


RW 


Rp/(Rp-X2) 


2  irk  z 

w  e 


in 


RW 


y<yX;> 


4.3.3  Pres£ure_drop  equ_i  1_i  bruim  i_n  the  heat  pipe 

Heat  transfer  in  a  heat  pipe  involves  the  working 
fluid  in  a  cyclical  process  in  which  the  fluid  is  pumped 
around  a  closed  loop  against  certain  pressure  losses.  During 
steady  state  operation,  there  must  exist  an  equilibrium 
between  the  pumping  forces  and  the  drag  forces  on  the  working 
fluid.  The  pressure  losses  undergone  by  the  fluid  arise  from 
the  viscous  resistance  to  the  flow  in  both  the  wick  and  the 
vapour  core,  and  from  the  momentum  loss  (or  gain)  of  the 
vapour  due  to  the  mass  transfer  at  the  liquid/vapour  interface. 


The  pumping  forces  are  the  capillary  force  on  the 
liquid,  which  is  due  to  the  variation  in  the  radius  of  curva¬ 
ture  of  the  liquid  menisci  along  the  wick,  and  the  body  forces 
on  the  fluid,  which  are  limited  to  the  force  of  gravity  in 
the  present  case.  The  weight  of  the  vapour  is  very  small  com¬ 
pared  to  the  weight  of  the  liquid  and  is  therefore  neglected. 


This  equilibrium  is  expressed  by  the  following  equa 
tion,  where  pumping  pressure  gradients  are  assumed  to  be  posi 
tive  in  the  flow  direction,  and  drag  pressure  gradients 
negative  : 
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APc+APg+aPe+APv  =  0 

where  aPc  pressure  difference  due  to  capillary  forces 

APg  pressure  difference  due  to  gravity 

ap^  pressure  drop  in  liquid 

ap..  pressure  drop  in  vapour 


(4.33) 


This  equation  can  be  written  in  integral  form,  with 
z  defined  as  the  direction  of  liquid  flow  : 


dpc 

% 

+  - 2-  + 

dp  ‘ 
V 

dZ 

dz 

dz 

dz 

dz  =  0 


(4.34) 


4.3.4  Pumping  forces 


In  the  case  of  a  class  B  wick  (in  which  the  reduced 
capillary  force  is  independent  of  the  wetting  angle  of  the 
liquid/solid  interface,  but  is  rather  a  function  of  pore  geo¬ 
metry),  the  maximum  capillary  pressure  gradient  is  given  Ly 


dz  r„(*«+£*  +  0 
p  e  a  c 


where  o  surface  tension  of  1 iquid/vapour  interface 

r  effective  pore  radius  in  the  wick 
P 


Gravitational  liquid  pressure  gradient 

This  term  is  simply  the  density  of  the  liquid 
times  the  acceleration  due  to  gravity  : 


4.3.5  Drag  forces 


Viscous  liquid  pressure  gradient 


This  term  is  given  by  a  form  of  Darcy's  law,  adopted 
for  flow  in  porous  media.  It  is  usually  used  for  fully  satu- 
tated  wicks,  with  A  denoting  the  full  cross-sectional  area  of 
the  wick.  In  this  case,  however,  A  represents  only  the  wetted 
area  and  depends  on  the  thickness  of  the  liquid  layer  in  the 
wick.  Cotter  (Ref.  4)  gives  the  following  expression  for  the 
liquid  pressure  gradient  : 

dpe  Mz) 

dz  K  A ( z ) 


where  m^ 

1  ocal 

liquid  mass  flow  rate 

(m/s) 

K 

wick 

permeabi 1 i ty 

(m2) 

A 

local 

liquid  flow  area 

(n>2) 

Vapour  pressure  gradient 

The  vapour  pressure  gradient  can  be  divided  into 
two  distinct  terms,  one  due  to  viscous  shear  and  the  other 
due  to  mass  transfer,  indicated  by 


'  dPv) 

’dp  1 

— -1  and 

V 

.dz  If 

.dz  . 

respectively.  The  expression  for  the  viscous  pressure  loss 
is  similar  to  that  for  the  liquid.  The  momentum  pressure 
gradient  depends  on  the  rate  of  mass  addition  (or  removal)  at 
the  wick  boundary  and  thus  contains  the  derivative  of  the  mass 
transfer  rate  with  respect  to  position.  Katsoff's  (Ref.  5) 
equations  for  each  of  these  terms  are 
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vvm(z) 


dP/ 

dz 


in 


8  m ( z )  dm (z ) 

3tt2pvRv  dz 


where  Ry  radius  of  vapour  core 

Vy  kinematic  viscosity  of  vapour 

Pv  density  of  vapour 


4.3.6  The  he^t_p/pe  equ_ations^ 


The  liquid  profile  in  the  wick  chosen  for  this  model 
imposes  certain  boundary  conditions  on  the  pressure  balance 
equations.  Since  the  thickness  of  the  liquid  layer  is  constant 
over  the  whole  condenser  section  and  the  whole  evaporator 
section,  according  to  the  resistance  model  the  heat  transfer 
rate  also  must  be  constant.  This  implies  that  in  the  pressure 
balance  equations,  ^  is  constant  over  each  section.  Additio¬ 
nally,  since  there  is  no  mass  interchange  with  the  reservoir, 
the  mass  flow  rate  is  zero  at  both  ends  of  the  heat  pipe. 


Expanding  equation  4.32,  and  integrating  4.34  over 
the  length  of  the  evaporator  nd  v.  /  the  length  of  the  con¬ 
denser,  one  obtains  the  following  >f  three  non-linear 

algebraic  equations  in  the  primary  and  secondary  unknowns 
introduced  in  section  4.21  : 


Mxi  >x2,x3)  =  0 
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(4.37) 

A  complementary  set  of  relations,  based  on  the  thermal 
resistance  model  resistance  model  of  figure  21,  is  used  to 
evaluate  the  secondary  unknowns  : 


Tv  =  Te  -  x3L (RPe+RWe) 


(Ta)  =  Tv  +  x 3L (RW  ) 


<VC  =  Tv 


x3L(RWc) 


(4.38) 


(4.39) 
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2 


(4.40) 
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4.3.7  Method_crf  s_o]_u_t i on 

In  the  FORTRAN  program,  equations  4.35  to  4.37  are 
solved  iteratively  using  Newton's  method.  The  input  to  the 
heat  pipe  subroutine  consists  of  the  known  heat  pipe  proper¬ 
ties  and  the  external  temperatures  of  the  evaporator  and 
condenser  sections.  The  output  consists  of  the  heat  flux  and 
heat  pipe  operating  temperature. 

A  matrix  is  constructed  of  the  partial  derivatives 
of  the  three  functions  fls  f2  and  f3,  defining  a  system  of 
linear  independent  equations  which  are  then  solved  by  Gaussian 
elimination  to  yield  a  correction  vector  containing  the  in¬ 
crements  to  be  added  to  each  of  the  primary  unknowns  in  order 
to  improve  the  accuracy  of  the  initial  estimates.  After  each 
iteration,  the  secondary  unknowns  are  reevaluated  using 
equations  4.38  to  4.40,  and  the  new  values  are  used  in  deter¬ 
mining  the  values  of  the  fluid  properties  to  be  used  in  the 
next  iteration.  This  procedure  is  repeated  until  the  desired 
tolerances  on  the  x  values  have  been  met. 

The  stability  of  this  method  seems  to  be  very  sen¬ 
sitive  to  several  of  the  physical  heat  pipe  parameters,  notably 
the  permeability  and  porosity  of  the  wick.  A  damping  factor 
is  currently  used  to  retard  convergence  and  thereby  to  avoid 
accidental  divergence  from  an  actual  solution. 

4 .4  Coupling  the  finned  tube  model 
to  the  heat  pipe  model 

The  critical  step  in  modelling  the  entire  heat 
exchanger  is  the  incorporation  of  the  heat  pipe  model  into 
the  previously  described  solution  of  the  finned  tube  system. 
What  follows  is  a  brief  explanation  of  the  global  solution 
algorithm. 


To  start,  an  overall  heat  exchanger  effectiveness  of 
50%  is  assumed,  and  the  heat  fluxes  in  all  of  the  heat  pipes 
are  assumed  equal.  This  gives  initial  values  for  the  heat 
fluxes  ^  in  vector  C  of  the  system  of  equations  shown 

in  figure  19.  The  inlet  air  temperatures  are  given,  so  the 
corresponding  set  of  equations  for  both  sides  of  the  heat 
exchanger  can  be  solved,  thereby  giving  the  air  and  fin  tempera 
tures  at  each  node.  These  temperatures  are  computed  on  the 
basis  of  an  arbitrarily  assumed  heat  flux  which  most  probably 
does  not  correspond  to  the  actual  performance  of  the  heat  pipes 
at  the  specified  air  inlet  conditions.  When  the  heat  pipe 
subroutine  is  called,  however,  new  values  are  provided  for  the 
heat  flux  in  each  row  of  heat  pipes.  These  new  values  replace 
the  previous  ones  in  the  vector  C_  in  the  parallel  systems  of 
equations  (evaporator  and  condenser  side),  enabling  one  to 
recompute  the  air  and  fin  temperatures  at  each  node.  The  pro¬ 
cess  is  repeated  until  convergence  is  obtained. 


The  output  values  of  particular  interest  are  T  for 
each  side  of  the  heat  exchanger,  which  correspond  to  TE(2)  and 
TC  ( 2 )  in  the  terminology  of  the  experiments.  With  these  data, 
the  heat  exchanger  effectiveness  is  computed  to  facilitate 
comparison  with  the  experimental  results.  The  operating  tem¬ 
perature  of  each  row  of  heat  pipes,  passively  determined  in 
the  process  of  convergence,  is  included  in  the  program's  out¬ 
put,  as  it  is  of  interest  in  determining  the  operating  regime 
of  the  heat  exchanger  predicted  by  the  program. 

4.5  Computational  results 


Since  we  had  almost  no  information  on  the  internal 
construction  of  the  heat  pipes,  informed  guesses  had  to  be  made 
to  provide  the  computer  program  with  some  of  the  necessary  in¬ 
put  parameters,  particularly  those  concerning  the  wick. 

Although  perhaps  not  the  real  configuration,  the  wick  was  assumed 
to  be  of  single  annular  geometry  with  isotropic  and  uniform 
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porosity.  Several  physical  properties  seem  to  play  a  major 
role  in  determining  the  stability  of  the  procedure,  notably 
the  porosity  (void  fraction)  and  permeability  of  the  wick. 

The  values  of  these  parameters  were  chosen  so  as  to  obtain 
convergence  of  the  solution;  they  do  lie,  however,  within!  the 
range  of  physically  reasonable  values.  A  summary  of  the  heat 
exchanger  parameters  used  in  the  computations  appears  in 
Appendix  C. 

The  program  HPIPE  was  run  for  a  variety  of  air  flow 
conditions  matching  those  of  the  experiments.  A  sample  of 
the  program  output,  with  only  the  final  iteration  shown, 
appears  in  Appendix  D.  The  air  and  fin  temperature  is  given 
for  each  node  (fin-tube  intersection)  as  well  as  the  vapour 
temperature  (reservoir  temperature  for  each  row  of  heat  pipes. 

A  plotting  routine  was  developed  for  the  program,  permitting 
the  visualization  of  the  evaporator  and  condenser  section  air 
and  fin  temperature  profiles  given  by  the  numerical  solution. 
Figures  22-29  show  the  temperature  profiles  computed  for  the 
various  air  inlet  temperatures  and  mass  flow  ratios  tested. 

It  is  important  to  note  the  difference  in  nature  between  the 
air  and  fin  temperature  profiles,  whereas  the  air's  tempera¬ 
ture  changes  monotoni cal ly  as  it  passes  through  the  heat 
exchanger,  the  fin  temperature  gradient  often  changes  sign  at 
the  heat  pipe  rows,  indicating  a  heat  flux  either  into 
(evaporator)  or  out  of  (condenser)  the  heat  pipes  both  upstream 
and  downstream  of  each  row.  The  existence  of  these  cusps  in 
the  fin  temperature  profile  depends  on  the  relative  magnitude 
of  the  equivalent  heat  pipe  thermal  conductivity,  the  fin 
thermal  conductivity,  and  the  coefficient  of  convection  of  the 
air  flows. 


For  the  highest  thermal  load  tested,  correspond!' ng  to 
an  inlet  temperature  of  65°C,  the  four  effectivenesses  computed 
for  the  four  different  mass  flow  ratios  are  in  fairly  good 
agreement  with  the  experimental  curves  when  plotted  versus 
their  respective  average  normalized  reservoir  (Fig.  30). 


Thus,  it  seems  that  despite  the  large  uncertainty  associated 
with  many  of  the  heat  pipe  parameters,  the  numerical  heat 
exchanger  model  does  approximate  the  heat  exchanger  effective¬ 
ness  in  one  particular  operating  regime.  There  is  a  major 
discrepancy,  however,  in  that  the  reservoir  temperature  pre¬ 
dicted  by  the  program,  which  is  based  on  the  assumption  of 
adiabatic  reservoirs,  are  consistently  higher  than  the  experi¬ 
mentally  determined  adiabatic  operating  temperature  for  the 
same  thermal  load.  This  discrepancy  can  perhaps  be  corrected 
with  a  more  detailed  knowledge  of  the  real  heat  pipe  parameters 
It  is  possible,  however,  that  the  assumption  of  a  constant 
liquid  profile  in  the  wick  of  the  evaporator  and  condenser 
sections  is  an  oversimplification,  inadequate  to  simulate  the 
actual  heat  pipe  behaviour,  which  may  include  significant 
dryout  zones  or  non-linear  liquid  thickenss  gradients. 

To  permit  the  numerical  simulation  of  the  heat 
exchanger  over  the  entire  normalized  temperature  operating 
range,  as  well  as  to  improve  the  prediction  of  the  adiabatic 
operating  point,  two  modifications  of  HPIPE  are  currently 
being  undertaken  : 

(1)  The  heat  pipe  simulation  will  be  expanded  to  take  into 
account  an  arbitrary  heat  input  to  the  reservoir,  hence  a 
vapour  flux  from  the  reservoir  to  the  heat  pipe.  For  given 
air  inlet  conditions,  the  operating  temperature  and  effective¬ 
ness  of  the  heat  exchanger  will  be  functions  of  the  reservoir 
input  energy. 

(2)  Rather  than  treating  the  evaporator  and  condenser  sections 
as  single  entities  with  constant  liquid  profiles,  each  section 
will  be  divided  into  finite  elements,  with  the  possibility  of 

a  different  liquid  thickness,  heat  transfer  rate,  and  tempera¬ 
ture  boundary  conditions  for  each  element. 
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5.  CONCLUSION 

The  research  presented  in  this  report  confirms  that 
the  use  of  heat  pipes  permits  the  construction  of  compact  and 
effective  air-to-air  heat  exchangers.  The  liquid  reservoir 
VCHT  can  be  considered  a  sort  of  thermal  transistor  (evapora¬ 
tor  =  emitter,  condenser  =  collector,  reservoir  =  base)  in 
which  the  primary  heat  flux  between  the  evaporator  and  conden¬ 
ser  is  modulated  by  a  secondary  heat  flux,  the  energy  input 
to  the  reservoir.  The  normalized  reservoir  temperature,  the 
mass  flow  ratio,  and  the  effectiveness  are  three  non-dimensio- 
nalizations  which  allow  the  coherent  representation  of  heat 
exchanger  performance  over  a  wide  range  of  operating  conditions. 
The  effectiveness  data,  when  normalized  in  this  way,  lie  on 
curves  which  exhibit  two  important  properties  : 

(1)  There  exists  a  peak  in  the  effectiveness  curves  which 
indicates  both  the  maximum  attainable  heat  flux  and  the  opera¬ 
ting  temperature  at  which  this  maximum  is  reached. 

(2)  Below  this  temperature,  there  is  a  nearly  linear  variation 
of  effectiveness  over  a  broad  temperature  band,  which  would 
facilitate  the  design  of  a  control  system. 

For  control  pupposes,  the  heat  exchanger  can  be 
considered  either  reservoir  temperatu re  modulated  or  reservoir 
power  modulated.  The  power  modulated  concept  is  more  direct, 
since  temperature  regulation  is  ultimately  achieved  by  power 
modulation  anyway. 

The  HPIPE  program,  once  it  has  been  modified  to 
simulate  complete  variable  conductance  operation,  will  be  a 
useful  tool  in  the  design  of  VCHP  heat  exchangers  and  in  the 
development  of  a  suitable  control  system.  It  appears  that 
the  structure  of  the  program  is  both  sound  and  efficient  and 
that  with  the  aid  of  more  complete  design  information  and  a 
discretization  of  the  heat  exchanger  along  the  length  of  the 
heat  pipes  as  well  as  in  the  flow  direction,  a  high  degree  of 
accuracy  can  be  obtained. 


ruiT^-THUrfU  I »  nm— C»| 
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APPENDIX  A  -  HEAT  EXCHANGER  SPECIFICATIONS 


Xp 

number  of  heat  pipes  : 

16  (in  a  square  4  > 

4  array) 

number  of  reservoirs  : 

4  (1  per  row  of  heat  pipes) 

Materials 

housing  •  aluminum 

fins 

:  aluminum 

heat 

pipes  :  copper 

reservoirs  :  reservoirs 

wick 

;  ? 

working  fluid  :  methyl ch 1 ori de 

(Freon  40) 

Dimensions 

evaporator  length 

.24  m 

e 

ir 

condenser  length 

.22  m 

c 

iU 

length  of  adiabatic  zone 

.02  m 

a 

S 

fin  spacing 

1.72  mm 

w 

width  of  heat  exchanger 

.  13  m 

depth  of  heat  exchanger  (flow  path  length) 

.13  m 

d 

diameter  of  heat  pipes 

.013  m 

P 

heat  pipe  spacing 

.  034  m 

diameter  of  reservoirs 

.020  m 
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APPENDIX  B 


SAMPLE  DATA  ACQUISITION  SYSTEM  OUTPUT 


DATE-' 7/5/31 
TITLE: .8/6/65 

CQNF I GURAT I  OH  COUNTER*  1 ,  PflRRLLEL=0  •'  I 


CAR ACT ER I ST I QUES  DES 

CANNAUX  A/D 

CANAL  1 

THERM 

CANAL  2 

ME  VAR 

CANAL  3 

MCOND 

CANAL  4 

NON  RACCORDE 

CANAL  5 

NON  RACCORDE 

CANAL  6 

NON  RACCORDE 

CANAL  7 

NON  RACCORDE 

CANAL  S 

NON  RACCORDE 

0 

TEMPS 

;:  0  SEC 

TEMP 

0-7 

+  .0 

+32. 2 

+32. 0 

+72.2  +72.8 

+67.  2 

+  .0 

+ 

0 

TEMP 

8-15 

+58.3 

+5y  a  y 

+58.8 

+58.2  +59.2 

+59.4 

+58.4 

+59. 

9 

TEMP 

16-23 

+52.  3 

+53. 1 

+52 . 8 

+52.3  +51.3 

+50 .  7 

+50.  7 

+58. 

6 

TEMP 

24-31 

+  .0 

+  .0 

+  .9 

+  .  0  +  .  0 

+  .0 

+  .8 

+ 

0 

TEMP 

•  "• !  -!  4 

+  .  0 

+  .9 

+  .  0 

+  .0  +  .0 

+  .  0 

+  .  0 

+  a 

0 

TEMP 

40-47 

+59.  1 

+59. 1 

+60.O 

+60.1  +58.9 

+61 . 8 

+68.8 

+  . 

0 

TEMP 

48-55 

+50.4 

+59 . 8 

+49.  8 

+58.1  +51.6 

+51.1 

+50. 8 

+ 

0 

PE: 

O  a  b  'Z'i*' 

696  MM 

1  H20 

MEVAP 

. 1 86228882 

:  KG/: 

SEC 

PC : 

6. 448 

71336 

MM  H2( 

It  MCOND  :  .  079252S 

:5Q1  1 

cG/S'EC 

TRES : 

67. 

2232388 

:  "C 

MRS: 

S  FLOW  RATIO’ 

. 746052747 

NORMAL I ZED  RESERVOIR  TEMPERATURE :  . 86878022 

TEG  >  •'  72.528554  TEC 2):  53.3699916 

T C C 1 >  :  32 . 1 439782  T C < 2  > :  51. 232 1 6 1 8 

HEAT  TRANSFER  RATE:  QE=  1.409O0195  KW 
POWER  RATIO:  QE/QC=  .921199642 
EE  •'  .  436749419  EX  •'  .  47 4 1 45963 

NORMALIZED  RESERVOIR  HEAT  TRANSFER:  QR/QE= 
LONGITUDINAL  TEMPERATURE  DISTRIBUTION 


QC=  1 . 52967886  KW 


•  88564 1'- 08 1''  1  OR."  Ql  = 


TLC 

1  >= 

60. 2520775 

TLC 

2  )  = 

59. 7690912 

TLC 

3  >  = 

59.1913034 

TLC 

4  >  = 

58 . 9864947 

TLC 

5  >= 

58. 6503912 

TLC 

h  ')  = 

51.2399422 

TLC 

"7  *| 

50.230832 

TLC 

g  )  = 

50. 6328261 

TLC 

9  ';»  = 

51 . 5049626 

TLC 

10  )  = 

52. 792086 

.  0788993:; 
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TEMPS'  109  SEC 


TEMP 

0-7 

+ 

.  0 

+31.9 

+31.7 

+72.0 

T  1  .  O 

+66.  7 

+  .8 

+ 

TEMP 

8-15 

+57.9 

+58.4 

+58.4 

+57.7 

+58.9 

+59. 1 

+58. 0 

+: 

TEMP 

16_2J 

+52. 1 

+52.7 

+52.4 

+51.9 

+58 . 3 

+5G.4 

+58.  3 

+: 

TEMP 

24-31 

+ 

.0 

+ 

.  0 

+ 

.  8 

+  .  8 

+  .0 

+  .  8 

+  .0 

+ 

TEMP 

32-39 

+ 

.  0 

+ 

.  0 

+ 

.0 

+  .8 

+  .  0 

+  .  8 

+  .8 

+ 

TEMP 

40-47 

+53 . 7 

+53 . 6 

+59.  6 

+59 . 8 

+58.5 

+60.  5 

+60. 4 

+ 

TEMP 

48-55 

+58. 1 

+58.5 

+49.5 

+49. 8 

+51.3 

+58.7 

+50. 6 

+ 

PE' 

8.  283 

184 

MM 

H20 

MEVAP : 

.  1  00 

1727553 

KG/S 

EC 

PC:  6.72312732  MM  H2G  MCOND •'  .0316234496  KG, -'SEC 

TRES:  66.7894924  "C  MASS  FLOW  RATIO :  .31033849 


SER' 

VOIR 

TEMPERATURE :  .  8658626; 

•025: 

i“r  •“» 

O  i  0 

TE<2> 

477 

OQ  i 
cLOH 

T  C  <  2  > 

50.' 9584425 

RATE : 

QE=  1 

.34238321  KW  1 

QE 

/QC= 

.  853 

502 

i 

■.1**1  *  ■ 

SER1 

VOIR 

HEAT 

TRANSFER •  QR/QE* 

EE:  .404252502  EX-  .473566619 

FORMALIZED  RESERVOIR  HEfiT  TRANSFER :  OF 
LONG I TUD I  NHL  TEMPERATURE  D I STR I BUT I  ON 


QC=  1 . 57245778  KM 


QR/QC5 


TLC 

1  >= 

59. 8800687 

TLC 

59.39658 

TLC 

3  y  - 

58. 7437604 

TLC 

4  >= 

58.6133442 

TLC 

5  >= 

53. 2776865 

TLC 

6  >  = 

58. 9365881 

TLC 

7  >= 

49.9868285 

TLC 

8  >= 

50. 3298574 

TLC 

9  •'= 

51.887948 

TLC 

10  >  = 

52.451878; 

TEMP 

'  8-7 

TEMPS'  2 
+  .8  +31 

TEMP 

3-15 

+53. 1 

+58.5 

+58.6 

+58.0 

+59. 1 

+3y  a  ^ 

+58.2 

+59.6 

TEMP 

16-23 

+52.2 

+52.9 

+52.7 

+52.0 

+51.1 

+56.  S 

+50.  4 

+50.7 

TEMP 

24-31 

+ 

.  0 

+ 

.  0 

+ 

.  8 

+  .8 

+  .0 

+ 

.0 

+  .0 

+ 

.0 

TEMP 

32-39 

+ 

.8 

+ 

.0 

+ 

.  O 

+  .8 

+  .  0 

+ 

.  0 

+  .  0 

+ 

.  0 

TEMP 

48-47 

+58.9 

+59.0 

+59.9 

+60 . 0 

+58. 9 

+68.  6 

+60. 5 

+ 

.0 

TEMP 

43-55 

+58 . 3 

+50.7 

+49 . 8 

+50. 1 

*51.5 

+51.1 

+50.  7 

+ 

.0 

PE:  8.6132592  MM  H20  MEVAP :  .10485466  KG/SEC 

PC:  6.72312732  MM  H20  MCOND:  .8316284496  KG/SEC 
TRES:  66.3492644  "C  MASS  FLOW  RATIO:  .784476633 


TEC  1 >  : 


SER 

VOIR 

TEMPERATURE :  .  85: 

:568! 

525 

TEC2> 

:  59.191121! 

10700 1 9 

TCC2> 

51.178915: 

;  RATE: 

QE=  1 

.  37842657  Kl 

QE 

/QC= 

.  364 

537943 

:454 

EX  ■  . 

477422173 

SER1 

VOIR 

HEAT 

TRANSFER :  Ql 

LONG I TUD I  NHL  TEMPERATURE  D I STR I BUT I ON 


TLC 

1 

>= 

60. 

866333 i 

TLC 

ii! 

)= 

59. 

6573855 

TLC 

>■' 

>= 

1*8 . 

967586 

TLC 

4 

■;i  = 

5b' . 

7^99645 

TLC 

5 

>  = 

58. 

464335 

TLC 

6 

)  = 

51. 

1641822 

TLC 

»’ 

'>  = 

50. 

2529397 

TLC 

U* 

..■>  = 

50. 

5189858 

TLC 

9 

>= 

51. 

2776361 

. 146366136 


.  0  +31.9  +31.8  •*■71.9  +72.5  +66.3  +  .0  +  .8 


QC=  1 . 58586325  KM 


i=  .15662834  QR/QC=  . 13541 2G58 
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APPENDIX  C  -  VALUES  OF  PARAMETERS  USED  IN  THE 
HPIPE  CALCULATIONS  HEAT  EXCHANGER 


Heat  exchanger 


length  of  evaporator  :  .24  m 
length  of  condenser  :  .22  m 
length  of  adiabatic  zone  :  .02  m 

Pipes 


external 

heat  pipe 

diameter  : 

.013 

m 

internal 

heat  pipe 

diameter  : 

.011 

m 

vapour  core  diameter  : 

.009 

m 

Wick 


wick  thickness 
poros i ty 
permeab i 1 i ty 

thermal  conductivity  of  wick  material 
effective  pore  radius 


.001  m 
15S 

1.  x  10-9  m2 
100  W/m°C 
.0005  m 


Fins 


fin  spacing 

.00172  m 

fin  thickness 

.0005  m 

thermal  conductivity 
length  of  platelets 

384  W/m°C 

.014  m 

X 

ro 

X 

u? 

X 

.034  in 

.014  m 
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Ut 


Working  fluid 


liquid  density  : 


Pn(T  )  =  -(1. 81 8x10~4)t]+ (1.125  10'2)T  -  ( 2 . 128  )T  +960 


£  '  V 


•V  ■  °c 


vapour  density  : 


pv  (Tv)  =  (1.441x10"“)t’-(5.762*10'3)Tv2+(.2965)Tv+6.066 


l  TVJ  • 


1 iquid  viscosity  : 


1  .  5 

p1  =  .6931  T  exp 


331.91 


t  T£  7 


Ns/m2 


It£]  -  0 K 


vapour  viscosity  : 


uv  =  (15.92  10"7)  T°‘5/  (1+462. 1/TV) 


Ns/m2 


[T0]  =  0  K 


latent  heat  of  vaporization 


4.2  x  10  J/kg 


thermal  conductivity  of  liquid  : 


.07 


W/m°C 


( kg/m  3) 


( kg/m  3 ) 
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APPENDIX 


-  SAMPLE  OUTPUT 


IP  I  PE 


*■»#***»•»»******»***»  *****»*»**»-s  #»*»*#  »*»*•■»•*#  »»#»4 •«•••• •*••••••••••••«•••••  •**•#««»*: 

PROPRIETE  DU  HEATPIPE 

RADP:.0055  KL:  .  07  KCU:  285.  RADV:  .  0045 

RADF:  0005  EPS:.  15  KM:  100.  MAXIT:50  RELHP:  .  2 

LE:  .  24  LC : . 22  LA:  . 02  RADO:  . 0065  KP: . 4648228692852 

L:  420000.  RPP:  .  00006905845561  127  Pl-JP  :.  34  22994  1 2G542 

Gi  1-5)  :  329450119.  7817  3.  183101  6504S9E+9  3.  1 0499 1  6054  1  5E+9  9.807  8333.333333333 
DELTA  ( l-3> :  1 .  E-10  l.E-10  l.E-10  N  3 

*•••••••«••««•»•••>•**«<•<«•<•••)«<<>••>•  *'»*»»*<.***-9  *•»*•»#»******«  «•»«**«*********«***»*»«*****■: 

FSO.  AIR  (A)  AND  FIN  <F>  TEMPERATURE  PROFILES  IN  THE 

EVAPORATOR  (E)  AND  CONDENSER  (C!  SECTIONS  OF  AN 
HEAT  PIPE  HEAT  EXCHANGER 
<+■)  HEAT  PIPE  VAPOR  TEMPERATURE 

AIR  INLET  TEMPERATURES  63  C,  30  C  ;  MASS  FLOW  RATIO  1. 

TEMP  (C) 

LENGTH  (M) 

TMAX  :  80.  TM  IN:  20,  NY:  6 
PROPRIETE  DU  CONDENSEUR 

MC:  .  75  TCI:  30.  HC :  42900.  LC:100.  F.C;.065  CPC;  1008. 

LCC: . 014  .  034  .  034  034  .  014 

PROPRIETE  DE  L 'EVAP03ATEUR 

ME.  75  TEI  65  HE:  46::00  LE  100.  EE’  .  07  CPE:  1008. 

LOE:  .  014  .034  034  .034  .014 

COUNTER:  T  MQwECH.  1 
XhP:  30. 

EPS:. 001  MLIM.  50  REL’  4  I MP I T  F 

##••##■*  ##«#####*##••*###«  *###».!#.?#*,*  »•••***#  *«##*•  **  #####*  #•##*##•*##«####*##■##«##..##♦*#  «»##*•#*' 

CONDENSEUR 

AC: 2. 049618074223  STC: 56  74603174603 

XLC :  7944444444444  1.729365079365  1  929365079365  1.929365079365  .7944444444444 
EVAPORATEUR 

AE:  1.  767272927273  STE’  61 .  11111111111 

XLE: .  8555555555556  2.  077777777778  2  0777777777’8  2  077777777770  .  0555555555556 


#  c 

*#•**•*#•*#••**#* 

r  -y  ♦  «  *  i  •*  • 

*****»*#*•#! 

» «  •  < 

** 

•*#######*#*•##• 

«*»«*« 

•*< 

»#•  < 

ITERATION  «0 

CONDENSEUR 

TA 

42  19 

41 .  24 

38. 

40 

35  40 

32. 

48 

31 

60 

T: 

42.  83 

43.  44 

4\) 

79 
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Mass 

flow 

ratio 

Heater 

power 

Normalized  reservoir 
temperatures 

12  3  4 

Ee 

ec 

1 . 

6  kW 

.37 

.33 

.36 

.32 

.29 

.30 

.8 

6 

.37 

.32 

.39 

.31 

.31 

.32 

.6 

6 

.37 

.32 

.39 

.30 

.37 

.36 

.4 

6 

.36 

.32 

.40 

.31 

.41 

.40 

1 . 

4 

.28 

.22 

.27 

.20 

.26 

.25 

.8 

4 

.32 

.26 

.32 

.26 

.30 

.30 

.6 

4 

.33 

.27 

.34 

.28 

.35 

.34 

.4 

4 

.31 

.26 

.33 

.28 

.42 

.37 

1. 

2 

.23 

.12 

.24. 

.14 

.22 

.21 

.8 

2 

.  19 

.  10 

.22 

.11 

.24 

.22 

.6 

2 

.16 

.07 

.18 

.09 

.27 

.19 

.4 

2 

.16 

.06 

.17 

.12 

.35 

.24 

TABLE  1 

-  SUMMARY  OF 

TEST 

RESULTS 

WITH 

INDIVIDUALLY 

HEATED  RESERVOIRS 

Row 

Mass  flow 

ratio 

T 

n 

c  e 

ec 

1 

.80 

.70 

.216 

.206 

2 

.80 

.70 

.218 

.224 

3 

.80 

.70 

.213 

.218 

4 

.80 

.70 

.185 

.216 

Heater  Power  =  4  kW 

TABLE  2  -  SUMMARY  OF  THE  TEST  RESULT  WITH 
INDIVIDUAL  RESERVOIR  BATHS 


variable  conductance  heat  pipes. 
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EMPERATURE  PROFILE  IN  THE  SU 


THE  HEA' 


CROSS 


/  Reserve Keoi.  (current)  source 
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THERMAL  RESISTANCE  NETWORK  MODEL  OF  THE  HEAT  PIPE 
(see  section  4.2.2) 


•IF:  A )  AMO  FIN  if)  TEMPERATURE  FROFILES  IN  THE 
EVAPORATOR  (E)  AND  CONDENSER  iC)  SECTIONS  OF  AN 
HEAT  PIPE  HEAT  EXCHANGER  . 


<  +  )  HEM  PIPE  VAPOR  TEMPERATURE 
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ATIONAL  AND  EXPERIMENTAL  RESULTS 


